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ABSTRACT 

We present a study of the dark and luminous matter in the isolated elliptical galaxy NGC 720, based 
on deep X-ray observations made with the Chandra and Suzaku observatories. The gas properties are 
reliably measured almost to R250O1 allowing us to place good constraints on the enclosed mass and 
baryon fraction (fb) within this radius (M25oo=l-6±0.2 x IO-'^^Mq, fb,25oo=0-10±0.01; systematic errors 
are typically <20%). The data indicate that the hot gas is close to hydrostatic, which is supported 
by good agreement with a kinematical analysis of the dwarf satellite galaxies. We confirm at high 
significance (^20-cr) the presence of a dark matter (DM) halo. Assuming an NEW DM profile, our 
physical model for the gas distribution enables us to obtain meaningful constraints at scales larger than 
R-2500, revealing that most of the baryons are in the hot gas. We find that fb within the virial radius 
is consistent with the Cosmological value, confirming theoretical predictions that a ~Milky Way-mass 
(Mvir=3.liJ5:3 X IQi^Mo) galaxy can sustain a massive, quasi-hydrostatic gas halo. While fb is higher 
than the cold (cool gas plus stars) baryon fraction typically measured in similar-mass spiral galaxies, 
both the gas fraction (fg) and fb in NGC 720 are consistent with an extrapolation of the trends with 
mass seen in massive galaxy groups and clusters. After correcting for fg, the entropy profile is close 
to the self-similar prediction of gravitational structure formation simulations, as observed in massive 
galaxy clusters. Finally, we find a strong heavy metal abundance gradient in the interstellar medium, 
qualitatively similar to those observed in massive galaxy groups. 

Subject headings: dark matter — Xrays: galaxies — galaxies: elliptical and lenticular, cD — galaxies: 
ISM — galaxies: formation — galaxies: individual (NGC720) 



1. INTRODUCTION 



In low-redshift disk galaxies of size comparable to the 
Milky Way, the fraction of cold baryons (cool gas and 
stars) relative to the total mass is typically found to be 
signi ficantlv less than the cosmological barvon fraction 
(e.g. iFukugita et al.||l998l iMcGaugh etldl bOlOl) , which 



stream (Moore fc Davis. 1994). gas stripping from dwarf 



satellites^ Blitz fc Robishaw 



2000) and provide a reser- 



is measured to be 0.17 (Dunkley et al. p009). This has 



voir for the ongoing gas accretion proposed to solve the 
"G-dwarf problem" (e.g. Sommer-Larsen et al.l^OOSl ). 

Direct observational evidence for such extended halos 
around disk galaxies remains controversial. While zero 
redshift absorption lines in the X-ray spectra of distant 



proven a challenge for standard models of galaxy forma- 
tion which, in fact, over-predict the observed barvon con - 
tent by as much as a factor ~2 (e.g. Benson et al. 2003 ). 
Solu tions to this problem, for example strong heating of 

the iiit.Prgtpllar mprliiiT ri (TC!\/n hy artivn cralar tir niirlpiis; 



the vicinity of the Milky Way (e. 


g. iNicastro et al. 


2002; 


Fanget al. 2002; 


Rasinusscii et al. 


2tX)3; b'ang et al. 


mi; 


Bregman & Lloyd-Uavies 


2UUV; Buote et al. 200(J|; 


bang 



(AGN) feedback (e. g. [Kauffmann et alj |1999| ) or less ef- 
ficient cooling (e.g. Mailer & Bullock 2004)7 ultimately 



line remains unclear (e.g. Yao fc Wan^ ^00"5| |Yao et al 
Other tracers of the putative hot baryons sug- 



2008) 



involi vo many of the baryonc remaining m a difficult to 



t Vint i p ttf 



mri P i t 



mrtnnrlnrl Vir 



detec 

around the galaxjf. ]\/Iodels suggest that as man}^ as 



gest the gas may be confined to t he disk ( Anderson 
Bregman 2010; Henley et al. 2010, although see loudaki 



l^ gjj et al. ' 2010 ), but the results are highly mo del- dependent 



the "missing baryons" in the local universe could be in 
diffuse, pressure-supported ha l os around normal galax- 
ies (e.g. I White fc Frcn"^ |1991|; iMaller fc Bul l(! ^ 12004 



Cerc p ot'' ^ft^ - 
,arsJ-r,|pnn(T| \ 



\ |b\ilaigita~fe Voobkuj ^j006 | ; | s 



Attempts to detect gaseous halos around external disk 
galaxies have, similarly, proven inconclusive. While X- 
ray emission from diffuse, hot gas associated with star- 
formation i n the disk has been o bserved (e.g. Fabbiano 



^iifli M half^ Miv^imrl ihn 1\TilVy Way nniilrl 



et al 



2001; Strickland et al. 2004), similar emission from 



help explain the prop erties of local OVI absorbing clouds 
( ^embach et al. 2003), the confinement of the Magellanic 



an extended, hot co r ona Has yet to be secu rely detected 



(Benson et al. 2000; Rasmussen et al. 2009). 
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In contrast to disk galaxies, dittuse X-ray emission from 
massive galaxy groups and clusters has been detected out 
to scales as large as ~Ri25o-R500^ > allowing tight con- 
straints to be placed on the total hot gas content (e.g 



Gastaldello et al. 2007b; Bun et al.| 2009; Vikhlinin et al 



^ We define as the geometrical radius within which the mean 
mass density of the system is A times the critical density of the 
Universe. 
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p006| ; [Allen et alj p008| ; [Pratt et al.|p009| ) . While the most 
massive clusters appear to be close to baryonic closure 
(i.e. the baryon fraction, fb, at large radii asymptotes to 
the to the Cosmologi cal value), fh at lower masses may be 
significantly smaller fpastaldello et aL 2007b| ; Pratt et al 
^009[ ; [Giodini et al.[ pDD9[ [Pai et alfl200g| ). The gas en 



tropy profiles of groups are systematically enhanced over 
the predictions of self-similar formation models, suggest- 
ing that non-gravitational processes (e.g. feedback during 
halo assembly) are responsible for e jecting barvons frorn 
the central regions of the halo fe.g. jPonman et al. 1999; 
[Pratt et alJ[2010HMcCarthy et al.[p010|) . Whether tEese 
baryons are completely ejected from the potential well 
is, however, uncle a r: extr apolating to Rvir, for example, 
[Gastaldello et al. (2007b) found fb close to the Cosmo- 
logical value for their group sample. 

Giant elliptical galaxies occupy an important niche be- 
tween Milky- Way sized disk galaxies and galaxy groups, 
offering a chance to connect these two mass regimes. As 
the likely end products of spiral galaxy merging, their 
baryon content may also provide important constraints 
on how the baryons are distributed in their progenitors. 
Unlike disk galaxies, most of their gas is typically in 
the form of an extended, hot halo, the X-ray emission 
from which is often detected out at least to tens of kpc. 
While early-type galaxies therefore provide an ideal op- 
portunity to take a census of the baryon content in nor- 
mal galaxies, surprisingly little is known about their fb. 
Based on weak lensing stu dies, the stellar mass fraction 
appe ars to be small (<0.03: Hoekstraet al. 2005; Gavazzi 
et al. ^007 ), but for massive galaxies, r nost of the baryon s 
are expected to be in the hot halo (e.g. Keres et al. 2005). 
In their hydrostatic X-ray mass analysis of th ree isolated 
ellipti cal gal axies (and four low- mass groups), Humphrey 
et al| ( 2006 , hereafter HOG ) were only able to obtain in- 
teresting constraints on the total baryon content by em- 
ploying restrictive priors . 

Mathews et aE (2005) showed that the X-ray lumi- 
nosity (Lx) of the most X-ray bright early- type galaxies 
is consistent with their being at the centre of massive, 
bary onically closed groups. However, there is observed 



an open question. 

To begin to address this question, in this paper we 
present a detailed X-ray study of the isolated elliptical 
galaxy NGC 720, using deep Chandra and Suzaku obser- 
vations (see Table |^) . The complementary characteristics 
of these satellites, i.e. the excellent spatial resolution of 
Chandra (enabling the inner parts of the X-ray halo to 
be studied in detail) and the low, stable background of 
Suzaku (which is helpful for studying the low surface- 
brightness outer parts of the halo), make a combination 
of these data ideal for studying the hot gas properties 
over as wide a radial range as possible. NGC 720 is a 
well-studied galaxy with a modest Lx; for its optical 
luminosity, it lies roug hly in the middle of th e scatter 
in the Lx-Lb relation (0 'Sullivan et al. ^001 ). It has 
occasionally been labelled a g alaxy group (e.g. [Garcia 
1993; Osmond & Ponman 2004), but within R500 the sys- 
tem comprises only one L* galaxy (NGC 720 itself) and 
dwarf companions, the brightest of whic h is three mag- 



nitudes fainter than the central galaxy (Dressier et al 
1986; Brough ct al. 2006). The total virial mass of the 



system, derived both from hydrostatic X-ray modelling 
and from the kinematics of the sa tellites, is a few times 
IO^^Mq ( pO^ ; [Brough ct al.[p006| ), making it only a few 
times more massive than the Milky Way, and therefore 
the ideal system for our purposes. Past studies, however, 
have not provided goo d cons traints on the overall baryon 
content of the system (H06). 

The X-ray emission from the hot gas in NGC 720, 
which can be traced out to scales of at least '-^BO kpc 




has a relax e d mo rphology ( Buote fc Canizares 



Buote et al 



a a huge range of Ly at fixed optical lummos- fc Valdarnini[[2008| ; [Fang et aL[p009[ ). The overaU agree- 
ity (L'anizarcs et ^ [19871 ; P'Sullivan et al.l t200l|; lEllisI ment between the masses inferred from hydrostatic meth- 
ods and those found by gravitational lensing and stellar 
dynamical modelling or the predictions of stellar popula- 
tion synthesis models generally support this p icture (e.g 



to bp a a huge range of Ly at fixed optical lummos 



2002), suggesting that hydrostatic 
mass modelling techniques can be reliably used to in - 
fer the total gravitating mass (e.g. [Buote fc Tsai[[l995( ). 
This is fu rther supported by th e lack of significant radio 
emission (Fabbiano et al, 1987), indicating that there is 
not currently an AGN that might be stirring up the ISM. 
For relaxed systems, hydrostatic equilibrium is believed 
to be an excellent approximation, with non-thermal ef- 
fects expected to contribut e no more than ^20 per cent 
of the total pressure (e.g 



fc Valdarnini 200? ; Fang et al 



>Jagai et 
009|) 



al. 2007; Piffaretti 



O 'Sullivan 2006 )7 suggesting that Lx is determined 
y both the virial mass of the halo in which galaxie s 
lie and their history of feedback (Mathews et al. 2006). 
In s ome cases, particularly in the lowest-mass systems. 



feedl ack appears to have becm strong enough to denude 



a galaxy largely of its gas ( pavid et al. p006 )^ Still, it 
is difficult to map galaxies from the optical versus X- 
ray luminosity plane onto fb directly, in large part since 
the gas emissivity depends on the square of its de nsit; 



Mahdiavi et alJ b008[ ; [Churazov ct al.[ ^008[ ; [Humphrey 
ct al. 2009a). Modest discrepancies have been reported 



odest discrepancies have been reported 
in a few galaxies between hyd rostatic X-ray masses and 



stellar dynamics results (e.g. Romanowsky et al. 200£; 



Gebhardt fc Thomai^p009| [Shen fc Gebhardt[[20iq), but 



and gas density profiles are not sclf-sirnilar (e.g. HO' 



nsity 
HOq: 



jCastjildello et al.( ^007b|; [Sun et all [20091). This meaSs 



these may simply renect systema tic uncertaint i es m the 
stellar dyna mical modelling (e.g. Gavazzi ^005 ; Thomas 



that Lx, which is typically measured within only a small 
fraction of the virial radius (Rvir), is a poor tracer of the 
overall gas mass. The detailed spatially resolved stud- 
ies needed to trace the gas distribution to large enough 
radii to measure fb more directly have largely been re- 
stricted to the most X-ray luminous systems, which tend 
to lie at the centres of massive groups. Thus, whether 
early-type galaxies hosted in Milky- Way sized halos can 
actually maintain a hot halo with a mass comparable to 
the predictions of disk galaxy formation models remains 



et al. 2007a ), or the dynamics of the hot gas at the very 
smallest scales (Brighenti et al. 2009] ). 

We adopted a nominal distance of 25.7 Mpc for 
NGC 720, based on the I-ban d SBF distance modulus 
estimate of Tonry et al. (2001), which was corrected to 
acco unt f or recent revisions to the Cepheid zero-point 
(see H06 ). At this distance, 1" corresponds to 123 pc. 
We assumed a flat cosmology with Hq — 70km s~^ and 
JIa = 0.7. We adopted R102 a s the virial radius f Rvir) , 
based on the approximation of Bryan fc Norman (l99S) 
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TABLE 1 

Observation summary 



ObsID 


Start Date 


Exposure (ks) 




Chandra 




7062 


2006 Oct 9 


23 


7372 


2006 Aug 6 


49 


8448 


2006 Oct 12 


8 


8449 


2006 Oct 12 


19 


800009010 


Suzaku 
2005 Dec 30 


177 



Note. — Details of the observations used in the present 
analysis. For each dataset we quote the observation identi- 
fication number (ObsID), the start date and the exposure 
time, after having removed periods of background "flaring" 
(Exposure). An additional, shallow Chandra observation 
(ObsID 492) was excluded due to its elevated background 
level. 

for the redshift of NGC 720. Unless otherwise stated, aU 
error-bars represent l-cr confidence hmits (which, for our 
Bayesian analysis, implies the marginalized region of pa- 
rameter space within which the integrated probability is 
68%). 

2. DATA ANALYSIS 

2.1. Chandra 

2.1.1. Data reduction 

The region of sky containing NGC 720 has been imaged 
by the Chandra ACIS instrument in the ACIS-S config- 
uration on 4 separate occasions (a shallow fifth observa- 
tion was excluded from our analysis due to its enhanced 
background level), as listed in Table J|. We processed 
each dataset independently, using the CTAO 4.1 and Hea- 
soft 6.8 software suites, in conjunction with the Chandra 
calibration database (Caldb) version 4.1.2. To ensure 
up-to-date calibration, all data were reprocessed from 
the "level 1" events files, following the standard Chan- 
dra data-reduction thread^. We applied the standard 
correction to take account of the time-dependent gain- 
drift and charge transfer inefficiency, as implemented in 
the CIAO tools. To identify periods of enhanced back- 
ground (which can seriously degrade the signal-to-noise, 
S/N) we accumulated background lightcurves for each 
dataset from low surface-brightness regions of the active 
chips, excluding obvious point-sources. Such periods of 
background "flaring" were identified by eye and excised. 
The final exposure times are listed in Table |l|. To com- 
bine the individual datasets, we merged the final data 
products (images, spectra, spe ctral response files, etc. ) 
using the procedure outlined in Humphrey et al] (2008). 
This involves correcting for relative astrometric errors by 
matching the positions of detected point sources (allow- 
ing for the possibility that some fraction of the sources 
are transient). Point sources were identified in each ob- 
servation by applying the CIAO wavdetect task to the 
full-resolution 0.5-7.0 keV band image, supplying the 
exposure-map (computed at an energy of 1.7 keV) to 
minimize spurious detections at the image boundaries. 
The detection threshold was set to 10~^, corresponding 
to < 1 spurious source detections per chip. A final source 
list was obtained by repeating the source detection pro- 
cedure on the final, merged image. All detected sources 



were confirmed by visual inspection, and, for each, appro- 
priate elliptical regions containing approximately 99% of 
its photons were genera ted. We detected 58 point s ources 
within the D25 ellipse ( |de Vaucouleurs et al. 1991), con- 
tributing ^^20% of the total 0.5-7.0 keV photons in that 
region. 

2.1.2. Chandra Image 

We examined the image for evidence of morphologi- 
cal disturbances that may complicate deprojection, or 
possibly suggest perturbation from hydrostatic equilib- 
rium. In Fig we show smoothed, fiat-fielded Chandra 
image, having re moved the point-s ources with the algo- 
rithm outlined in |Fang et al. ( [j009|) . The data were then 
flat-fielded with the exposure map and smoothed with 
a Gaussian kernel. Since the S/N of the image varies 
strongly with off-axis angle, the width of the Gaussian 
kernel was varied with distance according to an arbitrary 
power law, ranging from ^1" in the centre of the image 
to ~1' at the edge of the field. 

The image is smooth and slightly elliptical, showing 
no clear evidence of disturbances or asymmetries. This 
is consistent with our findings based on a shallower Chan- 
dra observation ( Buotc ct al. [2002|) . To search for more 
subtle structure in the image, we used dedicated soft- 
ware, built around the MINUIT library^, to fit an ellip- 
tical beta model (with constant ellipticity) to the central 
'^2.5' wide portion of the unsmoothed image, (correct- 
ing for exposure variations with the exposure map)^. We 
show in Fig ^ an indication of deviations from this sim- 
ple model by plotting (data-model)^ /model, which cor- 
responds to the residual in each pixel. To bring out 
structure, we smoothed this image with a Gaussian ker- 
nel of width 1.5" (3 pixels). As is clear from this "resid- 
uals significance" image, there is no evidence of any co- 
herent residual features. This is unsurprising given the 
lack of a currently active central black ho le, as suggested 



bv th e lack a significant radio emission (Fabbiano et al 
1987), nor a central, bright X-ray source. 



® http:/ /cxc. harvard.edu/ciao/threads/index. html 



2.1.3. Spectral analysis 

We extracted spectra in a series of concentric, con- 
tiguous annuli, placed at the X-ray centroid (which was 
computed iteratively and verified by visual inspection). 
The widths of the annuli were chosen so as to contain ap- 
proximately the same number of background-subtracted 
photons, and ensure there were sufficient photons to per- 
form useful spectral-fitting. The resulting annuli all had 
widths larger than '^8", which is sufficient to prevent the 
instrumental spatial resolution from producing strong 
mixing between the spectra in adjacent annuli. The 
data in the vicinity of any detected point source were ex- 
cluded, as were the data from the vicinity of chip gaps, 
where the instrumental response may be uncertain. We 
extracted products from all the active chips (excluding 
the S4 CCD, which suffers from considerable noise), mak- 
ing appropriate count-weighted spectral response ma- 
trices for each annulus with the standard CIAO tasks 

^ http: / /Icgapp.cern.ch/project / els /work- 

packages /mathlibs/ minuit / index.html 

* We obtained a best-fitting major-axis npro mrliiig r.f 9 R" j3 = 
0.38 and axis ratio of 0.85, consistent with Buote et al. (2002) 



Right ascension 



Fig. 1. — Top right: Suzaku XISO image in the 0.5-7.0 keV band, excluding data in the vicinity of the calibration sources. Circular 
exclusion regions (to mitigate bright point-source contamination) are overlaid. Bottom: Smoothed, point-source subtracted Chandra X-ray 
image in the 0.5-7.0 keV band. This image covers only the central part of the ACIS-S3 chip, where the count-rate is sufficiently high for 
the structure of the X-ray emission to be clearly discerned. The image contrast and smoothing scales were arbitrarily adjusted to bring 
out key features; the smoothing scale varies from ^l" in the central part of the images to ~1' in the outer regions. Logarithmically spaced 
contours are overlaid to guide the eye. Top left: "residual significance" image (see text) of the centre of the galaxy, indicating deviations 
from a smooth model fit to the X-ray isophotes. Although the isophotes are slightly elliptical, the galaxy appears relaxed and symmetric 
at all accessible scales given its distance. 



mkwarf and mkacisrmf . We extracted identical prod- 
ucts individually for each dataset, taking into account the 
astrometric offsets between them (determined from the 
image registration discussed above). The spectra were 
added using the standard Heasoft task mathpha, and the 
response matrices were averaged using the Heasoft tasks 
addrmf and addarf , with weights based on the number of 
photons in each spectrum. Representative spectra, with- 
out background subtraction, are shown in Fig 0, along 
with the best-fitting source and background models. 

Spectral- fitting was carried out in the energy-band 0.5- 
7.0 keV, using Xspec vers. f2.5.fn. We have shown pre- 
viously that fits to Poisson distributed data which mini- 
mize can yield significantly biased results, even if we 
rebin the data to more than the canonical ~20 counts per 

In contrast, we found that 



imizing C, but nonetheless took care to verify there is no 
significa nt bias, using th e Mont e Carlo procedure out- 
lined in Humphrey et al. ( 2009b ). Although not strictly 
necessary for a fit using the C-statistic, we rebinned the 
data to ensure a minimum of 20 photons per bin. This 
aids convergence by emphasizing differences between the 
model and data. The data in all annuli were fitted si- 
multaneously, to allow us to constrain the source and 
background components at the same time. 

We modelled the source emission in each annulus as 
coming from a single APEC plasma model with vari- 
able abundances, w hich was modified by Gala ctic fore- 
ground absorption ( Dickey fc Lockmaij 1990|) . Unlike 



several of our recent studies of galaxies ([Humphrey ct al 



bin ( Humphrey ct al 
the C-statistic of 



20091: 



"Fi~( r979 ) typically gives relatively jected spectra, similar to Gastaldello et al. (2007b|). We 



2006| , |2008i |2009at [Humphrey fc Buotc||20f 0| ) wc did not 



deproject the data, but instead fitted dirc c tlv the pro- 



unbiased results. We therefore performed the fits by min- 



adopted this procedure since we were interested in ac- 
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Fig. 2. — Representative Chandra and Suzaku XISO spectra for NGC 720, shown without background subtraction. In addition to the 
data, we show the best-fitting model, folded through the instrumental response (solid black line), along with the decomposition of this 
model into its various components. For the Chandra data, we show the hot gas contribution (solid red line), the composite emission from 
X-ray binaries (dash-dot magenta line), the instrumental background (dotted orange) and the cosmic X-ray background (dashed purple 
line). For the XISO data, wc show the hot gas emission from each annulus as solid lines (green from the central annulus, red for the second 
annulus and blue for the third), as well as the X-ray binary component, instrumental background and sky background (using the same 
colour scheme as for Chandra). For Chandra, the background is dominated by the instrumental component, whereas for Suzaku, which 
has a lower instrumental background but is less able to resolve the cosmic component into individual point sources, the cosmic component 
dominates. In all cases, emission from the ~0.5 keV gas is detectable above the background below ~1 keV. 



curately tracing the gas density to the largest available 
radii; deprojection requires the accurate subtraction of 
the emission from beyond the outermost radius, which is 
challenging given the very flat surface brightness profile 
of this galaxy. In our previous studies, our approximate 
method for achieving this introduced significant uncer- 
tainties into the outermost density data-point, which we 
have found to be prohibitive here. Since the temperature 
profile is relatively isothermal, the emission-weighted, 
projected gas density and temperature can be computed 
fairly accurately with the " response weighting " algorithm 
outlined in Appendix B of Gastaldcllo ct al.. We discuss 
the impact on our results ot fitting the deprojected data 
in S M. 



directly constrained. We allowed the abundance of Fe 
(Zpc), and the abundance ratios of O, Ne, Mg, Si and Ni 
with respect to Fe to vary. The other abundances (He) 
or abundance ratios with respect to Fe (for the other 



speci es) were fixed at their Solar values ( Asplund et al 
2004). To improve the constraints, Zpc was tied between 



To model the source spectrum, we adopted a slightly 



modified version of the Xspec vapec model, for which the 
abundance ratios of each species with respect to Fe are 



multiple annuli, where required, and the abundance ra- 
tios were tied between all annuli. The best-fitting abun- 
dances (Fig |) , are consistent (in the central part of the 
system) with previous Chandra results froin the shal- 
low (17 ks) observation ( Humphrey fc Buote| 2006| ), and 
are competitive with those obtained from the very deep 
Suzaku observation (§ 2.2; Tawara et al. 2008| ), as well as 
the XMM and Chandra measurements of the gas in much 
more massive ( and X-ray bright) objects (e.g. [Humphrcv 
fc Buote| p006| ; |Kim fc Fabbiano|p004| ; |Buote ct al.||2003| ) 



Intriguingiy, we note evidence of a significant abundance 
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Fig. 3. — Radial projected density, temperature and abundance profiles for NGC 720, obtained with botli Chandra and Suzaku. (For tlie 
definition of projected density, see text). Overlaid are the best hydrostatic models for each dataset, which match the data very well. The 
dotted lines are the best models if dark matter is omitted, for which the fit is very poor. The temperature and density data in the central 
bin of the Suzaku profile were excluded from the fit as the effects of the telescope's limited spatial resolution are most problematical at 
thiiLji£;ale, but are shown here as dashed diamonds. The central Suzaku abundance data-point appears sensitive to systematic uncertainties 
(§ Is.q) ; we show the results for a simultaneous fit of all the XIS units, and for a combination only of units 0, 2 and 3. The displayed 
error-bars do not take into account any covariance between data-points, which is accounted for in our analysis. In the upper right panel, 
we show the best-fitting, deprojected temperature profile, and the corresponding l-cr confidence range (shaded region). 



gradient in this system (discussed in more detail in § 3.5). 

To account for emission from undetected point sources, 
we included an additional 7.3 keV bremsstrahlung com- 
ponent. Since the number of X-ray point sources is 
approximately proportion al to the stellar light (e.g. 
Humphrey & Buotc ^008 ), the relative normalization of 
this component between each annulus was fixed to match 
the relative K-band luminosity in the matching regions, 
which we measured from the 2MASS image. To deter- 
mine the total X-ray binary luminosity, we added the 
total luminosity of this component, summed over all an- 
nuli, to the total Ly of the detected so urces within the 
D25 ellipse ( de Vaucouleurs et al. 1991 )^ The resulting 



Lx, (2.95 ±0.12) X 10 erg s~ , is in excellent agreement 
with that inferred from extrapolating the poin t source lu 
min ositv function ( ( 2 . 
& BT [otc||2008|) . 



± 0.8) X 10 "erg s : Humphrey 



To account for the background, we included additional 
spectral components in our fits, in a variant of the mod- 
elling procedure outlined in H06, which were fitted simul- 
taneously to each spectrum, along with the source. To 
account for the sky background, we included two (unab- 
sorbed) APEC components (kT=0.07 kcV and 0.2 keV) 
and an (absorbed) powerlaw component (F = 1.41). The 
normalization of each component within each annulus 
was assumed to scale with the extraction area, but the 
total normalizations were fitted freely. (We discuss the 
likely impa ct of the "Solar wind charge exchange" com- 
ponent in § 5.2.1). To account for the instrumental back- 



ground, we included a number of Gaussian lines and a 

^ To obtain Lx for the detected sources, we fitted an absorbed 
bremsstrahlung model to their composite spectrum. 



broken powerlaw model, which were not folded through 
the ARF. We included separate instrumental components 
for the front- and back-illuminated chips, and assumed 
that the normalization of each component scaled with the 
area of the extraction annulus which overlapped the ap- 
propriate chips. The normalization of each component, 
and the shape of the instrumental components, were al- 
lowed to fit freely. We included two Gaussian lines (at 
1.77 and ^2.2 keV), the intrinsic widths of which were 
fixed to zero. The energies of the ~2.2 keV lines were 
allowed to fit freely, as were as the normalizations of all 
the components. 

The best-fitting models are shown in Fig |^ for a rep- 
resentative selection of spectra. To verify the fit had 
not become trapped in a local minimum, we explored 
the local parameter space by stepping individual param- 
eters over a range centred around the best-fitting value 
(a nalogous to computing error-bars with the algorithm 
of pasl^p76| ). Er ror-bars were computed vi a the M onte 
Carlo technique outlined in Humphrey et al. ( 2006| ), and 
we carried out 150 error simulations. In our previous 
work these simulations were used to generate the stan- 
dard error on each measured temperature and density 
data-point, which were assumed to be uncorrelated. In 
a low-temperature (< 1 keV) system, however, there is 
a strong degeneracy between the abundance and the gas 
density. Tying the abundance between multiple annuli, 
as done in the present work, will consequently result in 
correlated density errors. Therefore, rather than com- 
puting only the error on each data-point from the error 
simulations, we also computed the covariance between 
each density data-point, whic h we fold into our mass 
modelling (see § ||). (In § 5^, we show that neglecting 
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this effect, or adopting a full treatment that incorporates 
covariance between all the temperature and density data- 
points, results in slightly different error-bars and a slight 
bias, but does not strongly influence our conclusions.) In 
Fig we show the derived gas temperature, abundance 
and projected density profiles and the l-u error-bars on 
each data-point (assuming the data-points are uncorre- 
lated). The "projected density" is the mean gas density 
if all the emission measured in the annulus originates 
from a region defined by the intersection of a cylindrical 
shell and a spherical shell, both of which have the inner 
and outer radii as the annulus. Defining it in this way 
is convenient as it has the units of density, but its pre- 
cise definition is unimportant, so long as care is taken to 
compute the models appropriately. 

As is clear from Fig H, the instrumental component 
dominates the background. This reflects both the strong 
non X-ray background of the telescope, but also the abil- 
ity of Chandra to mitigate the cosmic component by re- 
solving a significant fraction of it into individual point 
sources, which were excluded from subsequent analy- 
sis. Nevertheless, below ~1 keV, the emission from the 
~0.5 keV gas in NGC 720 dominates at least out to 
^60 kpc, allowing the gas temperature and density to 
be constrained. Emission from unresolved point-sources 
associated with NGC 720 is important only in the inner- 
most few bins, and (as is clear from Fig ||) it can be 
clearly disentangled from the hot gas component, as it 
has a very different spectral shape. We note that we have 
not explicitly included a component to account for the 
combi ned emission frorn cata clysmic variables and hot 
stars (Revnivtsev et al. 200S| ). This component, how- 
ever, is expected to be an order of magnitude fainter 
than the X-ray binary compon ent, and so this omissio n 
will not affect our conclusions ( Humphrey et al. 2009a ). 



2.2. Suzaku 

The observation of NGC 720 with Suzaku was made 
early in the life of the satellite, while all four XIS units 
were operating. Data-reduction was performed using the 
Heasoft 6.8 software suite, in conjunction with the XIS 
calibration database (Caldb) version 20090925. To en- 
sure up-to-date calibration, the unscreened data were re- 
pipelined with the aepipeline task and analysed follow- 
ing the standard data-reduction guidelinesfj. Since the 
data for each instrument were divided into differently 
telemetered events file formats, we converted the "5 x 5" 
formatted data into "3 x 3" format, and merged them 
with the "3 x 3" events files. The lightcurve of each in- 
strument was examined for periods of anomalously high 
background, but no significant amount of data was found 
to be contaminated in this way. In Fig |l|, we show the 
0.5-7.0 keV image for the XISO image, excluding data in 
the vicinity of the calibration sources. A number of off- 
axis point sources, most likely background AGNs, were 
identified by visual inspection. To minimize their con- 
tamination, we excluded data in a 1.5' radius circular 
region centred on each of these sources. Given the large 
telescope point-spread function, such a cut will only elim- 
inate ~70% of the photons from these sources. However, 
we have chosen not to use a larger exclude region, since 

^'^ http:/ /heasarc. gsfc.nasa.gov/docs/suzaku/analysis/abc/ 



it will cut out a large portion of the field of view while 
only slightly reducing the contamination. 

Spectra were extracted in three concentric annuli (0- 
2', 2-6' and 6-10'), centred at the nominal position of 
NGC 720 in the field of view of each instrument. Data 
in the vicinity of the calibration sources and the identi- 
fied point sources were excluded. For each spectrum, 
we generated an associated redistribution matrix file 
(RMF) using the xisrmfgen tool and an estimate of 
the instrumental background with the xisnxbgen task. 
Since the point-spread function of the telescope is very 
large (^2' half power diameter), even with such large 
apertures it is necessary to account for spectral mix- 
ing between each annulus. We did this by generat- 
ing multiple ancillary response files (A RFs) for each 



spect rum with the xissimarf gen tool (Ishisaki et al 



2007), which models the telescope's optics through ray 
tracing. Formally, we wish to model the true spec- 
trum of the hot gas, T^°'*(a, (5, £^), which is a function 
of right ascension (a), declination {5) and energy (E). 
To make the problem tractable, it is usual practice to 
adopt an approximation of the form T3'^^{a,S,E) ~ 
E,ar'(a,'5)ir'(-E)?»""(a,^), where af'(a,(5) simply 
delineates discrete, non-overlapping regions of the sky 
(i.e. a?"* (a, 5) = 1 if {a, 5) is within region i, or oth- 
erwise), and we have assumed that, within region i, the 
spectrum can be written tf"'*(£')^^°^(a, S), where tf°''^{E) 
is the normalized spectrum, and is the normalization 
as a function of position on the sky. This is analogous to 
expanding T^°-''{a,S, E) on a set of basis functions that 
are separable in E and (a, 5). For our purposes, each i 
refers to an annular region on the sky (with outer radii 2', 
6' and 10'), to correspond to the extraction regions. We 
approximated ^f"** as a /3-model (with parameters cho- 
sen to fit the Chandra image), correctly normalized to 
account for the gas emissivity. After passing through the 
telescope optics, a photon will be detected by the CCDs 
at some detector coordinate (x,y) and pulse height (h; 
corresponding to the energy). Th e pulse heigh t spectrum 
P{x,y, h) can be written as (e.g. Davis ^001 ): 



P(x, v,h)^ j da dS dEK{x, y, h, a, S, E)T{a, S, E) 

where K encapsulates the effects of the telescope optics 
and the instrumental response. Now, we define Dj{h) as 
the pulse height spectrum accumulated in an extraction 
region, j, such that Dj(h) = J dxdybj{x,y)P{x,y,h). 
(Typically we would choose i and j to correspond to the 
same set of regions on the sky, but that is not essen- 
tial). Here bj{x, y) is 1 if (x, y) lies within region j, or 
otherwise. Thus, we can write: 



dE 



J dx dy da d6 ar%j^''""'K 



Now, we assume that the term in square brackets can 
be replaced by a product of two new terms, Rj{E,h) 
and A^"*(£^). With this definition, these new terms cor- 
respond, respectively, to the redistribution matrix and 
ancillary response matrix that can be created with the 
tasks xisrmfgen and xissimarf gen. Hence, 

DAh) = E/^^ triE)R,{E,h)Af;^iE) (1) 
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Spectral fitting involves adopting a parameterized form 
for the spectrum ti{E) and comparing the corresponding 
model Dj (h) to the observed pulse-height spectrum. Un- 
like Chandra, for which the mixing between annuli was 
negligible (so that Afj'^{E) = if z 7^ j), it was neces- 
sary to to include all the ARFs Afj"* in the Suzaku fit- 
ting, which can be done with Xspec. The spectra from all 
extraction regions j were fitted simultaneously, to allow 
all three ti{E) models to be constrained. We note that 
the Xspec suzpsf model can be used to perform a more 
approximate correction for the mixing between annuli. 
However, our approach more formally accounts for the 
spatial distribution of the source photons, and is more 
practial for our present purposes. The spectra from all 
instruments were fitted simultaneously, allowing addi- 
tional multiplicative constants in the fit to enable the 
relative normalization of the model for each instrument 
to vary with respect to the XIS 0. 

It is necessary to add additional components to account 
for unresolved LMXBs and for the background compo- 
nents. The ARF terms {Aij{E)) depend not only on the 
extraction region geometry but also the spatial distribu- 
tion of the source photons, which clearly differs between 
the gas, LMXB and background components. Therefore, 
it is necessary to generate additional ARFs for each com- 
ponent. For each of these new components, we assume 
that the spectral shape does not vary with a and 6, and 
so we can write T'"''='\a, S, E) = ^''"^''(a, 
and similar terms for the background. We adopt a 
7.3 keV bremsstrahlung model for t^™^^[E), and allow its 
normalization to be fitted freely. Since LMXBs should 
trace the stellar light, we used the K-band 2MASS im- 
age (appropriately normalized) to define the term 
Folding all of the components into the fitting, eqn |^ is 
modified to: 

DAh)^Y. j tr{E)R,{E,h)A^;\E) 

i 

+ J dE t^"'''^{E)Rj{E,h)A]"'^\E) 
+ JdE t'^y{E)R,{E,h)Afy{E) 



+ / dE t"''\E)R,{E,h) 



(2) 



Here we consider separately the sky and instrumental 
("inst") background contributions. For the former, we 
adopted two AP EC pl asma models plus a powerlaw term. 



as described in § 2.1.3, and assume that its surface bright- 



ness does not vary over the field of view when com- 
puting A''^'^{E). For the latter term, we adopted two, 
complementary, approaches. One is to subtract off the 
spectrum generated for each annulus by xisnxbgen, in 
which case the last term in Eqn || will he omitted. The 
other method, which we adopted by default as it provides 
greater generality in our fits, is to model the instrumen- 
tal background, similarly to the treatment of this com- 
ponent with Chandra. We found that a single broken 
power law, plus Gaussian lines (at '^1.5, 1.7, 2.2, 5.9, 6.5 
and 7.5 keV), was sufficient for our purposes here, pro- 
vided the model was not multiplied by the effective area 
curve. Unfortunately, the quantum efficiency is folded 



into the RMF created by xisrmf gen, and so we first had 
to normalize the RMF through which this model was 
folded (the "renormalized" response matrix being indi- 
cated by Rj in eqn |^). Still, we found that the precise 
shape of the instrumental term is not important since 
it is a relatively minor component in the spectral model 
(except in the vicinity of the calibration lines). We did 
not obtain significantly different results if we used the 
modelling method or the direct sub traction method for 
the instrumental contribution (§ 5.2). 

We found that our model was able to capture the shape 
of the spectra well. The best-fitting normalization of the 
LMXB component corresponds to an LMXB luminos- 
ity of (3.08 ± 0.07) X lO'^'erg s~^, in excellent agreement 
with the Chandra data, and with the extrapolation of 
the detected LMXB luminosity function. Representative 
spectra are shown in Fig 0, showing the various model 
components. It is immediately clear that the spectral 
mixing between annuli is significant. For example, ^^20% 
of the photons from the central bin are scattered into the 
second annulus. Failure to take account of this effect ac- 
tually results in a significantly different temperature pro- 
file, which is almost isothermal at '^0.5 keV, and very 
inconsistent with the Chandra data. The background is 
dominated by the sky component, which is more signifi- 
cant than in Chandra since very few of the point sources 
that dominate it above ^1 keV can be resolved and ex- 
cluded. 

The hot gas temperature, abundance and projected 
density profiles are shown in Fig ^. In the two outer an- 
nuli, the data are in good agreement with the results from 
Chandra. In the innermost annulus, however, we found 
some discrepancies. Specifically, the abundance is sig- 
nificantly higher than observed with Chandra, while the 
temperature is slightly higher than an emission- weighted 
average of the Chandra temperature profile in this re- 
gion. The Suzaku results in the central bin are, how- 
ever, sensitive to systematic uncertainties in our mod- 
elling procedure (sec § |^). In particular, if we omitted 
data from the XISl unit from our fits, we obtained Zpo 
in much better agreement with Chandra (Fig H) . While 
this may point to problems with the calibration of the 
XISl unit, it may also be related to errors in the correc- 
tion we applied to account for spectral mixing between 
the different annuli. If we omit this correction, for ex- 
ample, the discrepancy with Chandra is significantly re- 
duced (Zfc~0.9). The real temperature profile is clearly 
more complex than we assume when performing this cor- 
rection, which may be partially responsible, as there may 
be slight errors in the ray-tracing algorithm. Further- 
more, the single-temperature model fitted to the large 
Suzaku regions, within which the Chandra data reveal 
temperature gradients, would not be expected to yield 
perfectly consistent results with the Chandra data. In 
any case, we omit the data from the central Suzaku annu- 
lus in our subsequent modelling, since the Chandra data 
are expected to be more reliable on these scales. 

3. MASS MODELLING 

We translated the density and temperature profiles 
into mass constraints using th e ent ropy -based "forward- 



fitting" technique outlined in HOB and Humphrey et al. 
( 2009a ) . Briefiy, given parametrized profiles of "entropy" 



Fig. 4. — Left: K-band 2MASS image of NGC 720. Centre: residual image, having subtracted our projected model for the three- 
dimensional light distribution. Overall, this model provides a good fit to the stellar light distribution. A small, coherent residual feature 
is visible towards the top of the image, which may not be associated with the galaxy. If it is a part of NGC 720, its luminosity is too faint 
to affect our conclusions. Right: the spherically averaged mass density as a function of 3-dimensional radius from our model. 

-2/3 



{S—kTue , where rig is the electron density) 
itating mass (excluding the gas mass, which is 
self-consistently) , plus the gas density at some 



radius (for which we used 10 kpc), the three-di nensional imag3 prior to deprojection. We show in Fig ^ the K 



temperature and density profiles can be calcula ,ed, under 
the hydrostatic equilibrium approximation and fitted to 
the data. This model assumes spherical symme ]ry, which 
is a standard assumption in X-ray hydrostatic i Qodelling, 
even when the X-ray isophotes are not perfec 
since deviations from sphericity are expecte 
tribute only a rel ativclv small error into the 



mass profile (e.g. Piffaretti ct al, 2003: Gava 



ly round, 
i to con- 
recovered 
20051), 



( 200! la), we subsequently ave raged it spherically tor us e 
with our X-ray modelling (see Buote & Humphrey 201C| ) . 

^ . To fit the three-dimensional entropy prohle we assumed 

and the inferred baryon fraction (Buote & Humphrey a model comprising a broken powerlaw, plus a constant; 
^010 ). such a model is reasonably successful at reproducing the 

~ ' ' entrc py profil es of galaxies and galaxy groups over a wide 



For the gravitating mass model (in addition bo the gas 
mass), we explored two possible models. Firstly^, we con- 
sider ed a model without dark matter, so that the grav- 



itatii Lg mass comprised only a stellar compon ant, a su- 
permassive black hole (with mass fixed at 3 k WMq, 
which is that estima ted from the Mbh-c* rslation of 
Tremaine et al. ^002 , given the central velocity disper- 
sion value listed in the HyperLEDA database , and the 
mass of the gas. We assumed the stellar mass c Dmponent 



followed the optical light, but allowed the ma 



(M/L) ratio to be fitted freelv. In the second nodel, we 
also included an NFW ([Navarro et alj |l997| ) (Hark mat- 
ter halo, the virial mass and concentration (defined as 
the ratio of the virial radius to the characteristic scale 
of the NFW model) of which were allowed to 
To model the stellar light, we deprojected th e K-band 
2MAS S imag e using the method outlined in 1 [umphrey 
et al [j 2009a ), which is based on the algorith n of Bin 
ney ('t aT ( p"990 ). This entailed first fitting an arbitrary, 
analytical model to the 2MASS image, which -sras depro 
jected analytically onto an arbitrary axisymmstric grid, 
for a given inclination angle. Since NGC 720 sxhibits a 
very elliptical optical morphology, it is improbable that 



it is being observed at a low (face-on) inclinat 



and so we initially assumed a 90° inclination ( j ut in ves 
tigate th e effe c t of u sing lower inclinations in § 5.6). A 



series of Lucy (1974) relaxa tion iterations gr a dually im- 
proved the deprojection. In Humphrey et al. (S009a), we 
demonstrated that this algorithm gives results consistent 
with the alternative deprojection technique of Gcbhardt 



and grav- 
omputed 
canonical 



et al 
extr; 
not I 



band 2MASS image, the residuals from the projection 
of or r axisymmetric model (demonstrating the excellent 
and the spherically averaged stellar light density 
c. Although this method yields an axisymmetric 



fit), 
profi 

three -dimensional light profile, following Humphrey et al 



radial range ( | tlumpnrcy ct al.| |z!UUija|; iJetna et al.||^UUY|; 
Fino uenov et al.| p007| ; [Gastaldcllo ct alj p007at [Sun 



et al 2001; The normalization of the powerlaw and 
cons' ant components, the radius of the break and the 
slope s above and below it were allowed to fit freely. We 
inves tigate the impact of varying the slope outside the 



s-to-light 



vma. 

onic 

viria. 

mass 

sible 



ion angle, dens: 



( 1996 ). Since the innermost region within which we 
cted the Chandra data has a radius of ^^10", it was 
ecessary to deconvolve the seeing from the 2MASS 



1 range (Humphrey ct al. 2009a ; Jetha et al. |2007 



; of the data (~80 kpc) between reasonable limits, 
1 shoul d bo und all plausible extrapolated entropy 
es, in S 



rangi 
whic. 
profi 

Fo lowing Humphrey et al. (200S), we solved the equa- 
tion I )f hydrostatic equilibrium to determine the gas prop- 
erties as a function of radius, from 10 pc to a large radius 
outside the field of view. For the latter, we adopted the 
radius of the system defined by ignoring the bary- 
;omponents (which is slightly smaller than the true 
radius of the system)]^ For any arbitrary set of 
and entropy model parameters, it is not always pos- 
to find a physical solution to this equation over the 
full radial range. Such models were therefore rejected as 
unpl ysical during parameter space exploration. In order 
to cc mpare to the observed projected density and tem- 
pera ure profiles, we projected the three-dimensional gas 
ty and temperature, using a procedure similar to 



br the constant M/L ratio model, we found it more conve- 
nient ;o use a less conservative, smaller radius. Specifically we used 
75 kp ;, and found the model agreed even less well if a larger scale 
was a iopted. 
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Fig. 5. — Best-fitting entropy profile from fitting the projected 
data (soUd fine). The grey shaded region indicates the approximate 
l-cr scatter. The sharpness of the break at ~14 kpc is a consequence 
of the simple, broken powerlaw parameterization we adopted for 
the entropy profile. Overlaid are data-points obtained directly 
from-the deprojected Chandra (triangles) and Suzaku (stars) data 
(§ |5.4[ ), which agree well with the projected fit. The deprojected 
data points have much larger errors as deprojection tends to am- 
plify noise, and the model fit uses both Chandra and Suzaku data. 



that described in pastaldello et alJ ( |2007bD p|. This in 



volves computing an emission- weighted projected mean 
temperature and density in each radial bin, that can be 
compared directly to the datapl 

To compare the model to the data, we used the 
statistic, fitting simultaneously the Chandra and Suzaku 
temperature and density profiles, with the central bin 
of the Suzaku data excluded from our fits. Correlations 
between the density errors were simply implemented by 
adopting a form for the statistic wh ich incorporates 
the covariance matrix (e.g. |Gould| |00|)[]- The proba- 
bility density function (of the data, given the model) is 
modified in a similar manner when folding in the non- 
diagonal elements of the covariance matrix. Parame- 
ter space was explored with a Bayesian Monte Carlo 
method. Spe cifically, we used version 2.7 of the Multi- 



Specu 

Nest cod«2l ( Feroz fc Hobson| ^008| ; |Feroz et alj ^0091 ). 



Since the choice of priors is nontrivial, we followed con- 
vention in cycling through a selection of different priors 
and assessing thei r impact on our results (we discuss this 
in detail in § 5^ ). Initially, however, we adopted fiat 
priors on the logarithm of the dark matter mass (over 
the range 10^^ < M < lO^'^M©), the logarithm of the 
dark matter halo concentration, cjjm (over the range 
1 < cdm < 100), the logarithm of the gas density at the 

In computing the plasma emissivity term, we approximated 
the true three dimensional abundance profile with the projected 
abundance profile (Fig since the profile is relatively flat, this 
should be sufficient for our present purposes. 

^'^ We note that, for systems with a strong temperature gradient 
and kT ^ 3 keV the emissi on- weighted temperature may be biased 

craiure of NGC 720 the 
iOOC ) , which we conflrm 
ITa when_fitting the de- 



20M). For tbc_tem 

DC weak (Buotc 



5.4) 



low (Vlazzotta et al. 
effect IS expected to 

given the lack of any significant bias we 
projected data, rather than the projected data ( 

By default we only consider correlations between the density 
data, but we invostigate introducing a more complete covariance 
computation in § |5.q and find that the results are not significantly 
affected. 

http:/ /www. mrao.cam.ac.uk/software/multinest/ 



TABLE 2 

Entropy and pressure fit results 



Parameter 



Marginalized value 



Best-fit 



units 



logio Ps.io 


-26.27±0.02 


(-26.28) 


[g cm 3] 


Sc 


3.9±0.6 


(3.7) 


keV cm^ 


So 


23±2 


(23) 


keV cm^ 




1.14±0.19 


(1.13) 




Rbrk 




(10.9) 


kpc 


/32 


0.47±0.13 


(0.51) 





Note. — Marginalized value and l-cr errors on the depro- 
jected gas density at 10 kpc (pg,io), and the entropy param- 
eters. The fitted entropy model within 100 kpc has the form 
S{R) = Sc + Sof{R), where f{R) = ijf^ , if _R < R^rk, or else 

f{R) = {Rbrk/iOkpcf'^-f^^Rl^, and i?io = i?/10 kpc. Since 
the best-fitting set of parameters need not be identical to the 
marginalized values, we also list the best-fitting value for each 
parameter. 




Fig. 6.— Radial mass profile of NGC 720. The solid (black) 
line indicates the total enclosed mass, the dashed (red) line indi- 
cates the stellar mass, the dotted (blue) line is the dark matter 
contribution and the dash-dot (magenta) line is the gas mass con- 
tribution. The grey shaded regions indicate the l-cr error on the 
total mass distribution. Overlaid are a set of data-points derived 

frorn the mnre "traditional" rriass analysis of this System described 

in Humphrey & Buotc (201C), which generally agrees very well 
witn tne fttted model (we stress the model is not fitted to these 
data-points, but is derived separately). Note that the kink in the 
stellar mass model at ~20 kpc is an artefact of the stellac, light 



deprojection technique, but does not affect our results (§ 5.6) 

canonical radius, the stellar M/Lk ratio, and the vari- 
ous entropy parameters. For a r nore detailed desc r iption 
of the modelling procedure, see Humphrey et al. ( 2008 , 
2009aD . 



4. RESULTS 

4.1. No dark matter model 

We first consider the case without a dark matter halo 
component. Although our previous work (e.g. H06| ) has 
found that the M/L ratio rises steeply with radius (im- 
plying that dark matter is needed at high significance), 
given the high-quality data in NGC 720, it is of partic- 
ular interest to quantify how poorly the constant M/L 
ratio model fits these data. We find that the best- 
fitting model without dark matter is a very poor de- 
scription of the temperature and density profiles (Fig ^; 
X^/dof=413/19, with a mean absolute fractional devia- 
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Fig. 7. — Radial total mass-to-K-band light profile of NGC 720. 
The shaded region indicates the l-tr scatter about the best-fitting 
relation. The steep rise with radius outside ~Rc indicates the pres- 
ence of a massive DM halo. The small kink in the profile around 
~20 kpc is an artefact of the stellar light deprojection technique, 
and is unrelated to the break in the entropy profile at ~10 kpc. 
Overlaid are a series of data-points derived from the more tradi- 
tional mass analysis shown in Figfel We stress that the data-points 
are derived independently and the model is not fitted to these data. 



tio n between the dat a and models of '^29%). As found 
by Buote et al. ( 2002 ), based on an analysis of the X-ray 
image in NGC 720, the stellar mass component is too cen- 
trally concentatcd to produce the observed gas pressure 
profile. 



4.2. Dark matter halo 

As is clear in Fig |[ when a dark matter halo com- 
ponent was added, the hydrostatic model was able to 
capture the overall shape of the projected density and 
temperature profiles very well (x^/dof=12.7/17). The 
improvement to the fit when dark matter was included 
was therefore highly significant; the ratio of the Ba yesian 
"evidence" retu rned by the fitting algorithm (see Fero^ 
[fe Hobson|[2008| ) for the two cases is 10 implying that 
dark matter is requiredFj at ~19.6-cr. 

In Fig 1^ we show the oest-fitting entropy profile (and 
its l-cr confidence range), which shows a flattening at 
large radius similar to t wo of the three galaxies stud- 



iedin Humphrey et al. ( [2009a ). As we will show later 
(§ 6.4), the entropy profile shape, even when extrapolated 
to large radius, is consiste nt with trends seen in galaxy 
clusters ( |Pratt et all|2010|) . The derived radial distribu- 
tion of the gravitating mass is shown in Fig o and the 
total mass-to-light ratio is shown in Fig ^ In Fig ^ we 
also show the contribution to the mass from each sep- 
arate (xmponm^^ 



mass distribution within the central ~5 kpc, but a signif- 
icant aSJFlnalleFTiatoliTeqTnred"^ 
at larger radii, consistent with our previous study of this 
system, using data from a much shallower observation 



^® Although not correct in this case (Protassov ct al 
reference, it is interestin] 



2002| ), for 
the dark 



to quote the significance o: 
matter detection in frequentist terms by using the F-test; in this 
case, such an improvement in would occur by random chance 
with a probability of 1.4xlO~^ , implying a ~7.4-a- detection of 
dark matter. 




12.4 12.6 
l°ilo(M™-/M0) 

Fig. 8. — Inferred Cvir-Mvir relation for NGC 720. We show the 
marginalized l-cr joint confidence region for Cvir and Mvir based 

(light lalue) curve is the 
L_and the dot-dash-dash 
lOq. The solid black line 



on our fits (solid blue line). The dott 
same but for the deprojected data (§ p.i 
region is the l-cr error region found byl 
and grey shaded 



scatter found by Buote et al 
the theoretical relation tor re 



rnginn arn 



o moan Cvir-Mvir 



relation, and l-cr 



(2007), and the .red shadrd rpmnn .is 
axed nalos from Maccio et al. (2008). 



(Humphrey et al. 2006). The hot gas becomes the dom- 
inant baryonic component outside '^200 kpc (^Rsoo)- 
Overlaid are a series of mass data-points derived from 
fitting the deprojected Chandra data using a more tradi- 
tional (but more unce rtai n) mass-modelling meth od (for 
more details see § 5^ and Humphrey et al. 2009a ), which 
agree well with the inferred mass distribution, indicating 
that the resulting mass profile is not overly sensitive to 
the analysis method. 

A notable feature of the recovered mass distribution is 
its fiattening outside ^20 kpc, which allows us to con- 
strain the scale radius, and hence the virial mass of the 
dark matter halo. In Fig ||, we show the relation between 
the concentration of the gravitating mass, Cvir, and the 
virial mass. M,,;,-. To be consistent with our past work 



( Buote ct al. 2007), the Mvh- and Rvir are derived from 



the distribution of the total gravitating mass, not just 
the dark matter, and the concentration, Cvir is defined 
as the ratio of Rvir to the characteristic scale of the DM 
halo. These results are consistent with those obtained 



by HOC with Chandra, but the constraints are very much 
tighter, refiecting both the significantly deeper observa- 
tion (~5 times deeper Chandra observation, combined 
with the data from the very deep Suzaku observation) as 
well as improvements in our mass modelling procedure. 
Overlaid on Fig |^ we show the theoretical concentration 
versus mass relation from dark matter only simulations 
dMaccio et"al| |2008|, adopting their "ACDMl" r elation 



for re l a.xed h alos ) , and the empirical relation from Buote 



et al 



( ^0071 ). NGC 720 lies abov e both relati ons, but is 
within the 2-cr scatter obtained by p3uote et~aL . The best- 
fitting marginalized model parameter values and confi- 
dence regions are given in Tables ^ and 0, 



In c ontrast to our measurement in Humphrey et al 
(2006), our adopted modelling procedure did not require 



us to adopt a strong prior on fb in order to obtain inter- 
esting virial mass constraints. This allows us to measure 
fb directly from our fit. In Fig ^, we show how the en- 
closed baryon fraction and gas fraction vary as a function 
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TABLE 3 
Mass results and error budget 



Test 


M,/Lk 


log M2500 


log C2500 


log M500 


log C500 


log Mvir 


log Cvir 




MqLq-i 


[Mq] 




[Mq] 




[Mq] 




Marginalized 


'-'■^^-0.05 




0.86 ±0.08 




1.15 ± 0.07 




1.43 ±0.08 


Best-fit 


(0.54) 


(12.20) 


(0.85) 


(12.36) 


(1.14) 


(12.51) 


(1.42) 



ADM 
AAC 

ABackground 
ASWCX 
AEntropy 
A3d 

ASpectral 
AWeighting 
AFit priors 
AInstrument 
AStars 
ADistance 
AFit radius 
ACovariance 



+0.07 (IH? 
-0.165 (±0.04) 
-0.077 (±0.05) 

-0.057 (±0.05) 

-0.014 (±0.04) 

-0.008 (±0.05) 

±0.04 (±0.06) 

-0.011 (IHO 
1^2 (±0.05) 
IH2 (±0.05) 
+°-^^ (±0.33) 



-0.03 



±0.04 
+0.01 
-0.019 (±0.06) 



-0.04 
-0.05 
-0.05 
-0.04 



+0.12 (±0.03) 
+0.03 (±0.06) 
-0.051 (±0.04) 

-0.059 (IH^) 
+0.005 (±0.06) 
-0.016 (±0.07) 

+0.17 ^'+0.08^ 
-0.03 1^-0.06 j 

+0.01 (±0.05) 
lao5 (±0.06) 
+0.10 

+0.00 /'+o.o8^ 

-0.01 1^-0.06 j 

iao? (±0.06) 
+0.09 (±0.11) 

+0.01 /'+o.o6^ 

-0.01 1-0.07 j 



+0.46 (±0.08) 
-0.092 (±0.09) 

+0.11 

+0.10 (±0.08) 
±0.01 (±0.08) 
+0.009 (±0.10) 
lo:o3 (±0.09) 



+0.06 



/+o.io^ 
U0.12J 



1^2 (±0.09) 
-0.087 (iHg 
+0.01 (iSi^ 
law (±0.08) 



error estimates 
+0.30 (±0.03) 
+0.05 (±0.07) 
-0.059 (±0.05) 

-0.055 (±0.05) 

+0.008 (±0.07) 

-0.031 
iaJs (±0.08) 
+0.01 (±0.06) 
tlZ (±0.07) 



+0.52 (±0.08) 
-0.087 (±0.09) 
+0.10 (±0.07) 

+0.10 (ino 

torn (±0.08) 
+0.01 (±0.10) 
(±0.09) 



hO.48 (±0.03) 
[-0.05 (±0.06) 
-0.059 



+0.05"! 
-0.04 1 
+0.05^ 
-0.03 1 
+0.091 
-0.06 I 
+0.08 
-0.04 



-0.11 



+0.11 



-0.072 



z'+o.oo^ 



V-0.127 
-0.015 (±0.10) 



-0.09 
+0.00 (+0.12 
-0.02 ^^-0.06 

tHt (±0.06) 



+0.10 (±0.12) 
+0.02 



fO.06 1 
-0.08 I 



-0.03 

+0.06 

+0.04 
-0.12 

-0.087 
+0.01 

io:o4 (±0.08) 
-0.068 rt»:f2 

-0.020 



+0.09 1 
-0.12 I 



h0.07 
-0.09 
+0.17 
-0.25 



-0.078 
-0.033 
-0.031 
taol (±0.08) 
+0.006 (±0.06) 
tHI (±0.06) 

+0.11 



+0.58 (±0.08) 
-0.089 (±0.09) 
+0.10 (±0.07) 

+0.10 (±0.08) 
1^2 (±0.09) 
+0.02 (IHI) 
lo:o4 (±0.09) 
+0.06 



0.018 {tin 

+0.06 
-0.05 
+0.13 
-0.10 



±0.03 
+0.09 



V-0.127 
(±0.08) 
-0.086 (±0.09) 

^+0.17^ 

(±0.08) 



(±0.06) 



-0.04 

-0.076 
-0.022 



'+0.09 1 
,-O.llJ 
'+0.121 
, -0.08 I 



Note. — Marginalized values and 1-a confidence regions for the stellar mass-to-light (M*/Lk) ratio and the enclosed mass and concentration measured 
at various ovcrdcnsities. Since the best-fitting parameters need not be identical to the marginalized values, we also list the best-fitting values for each 
parameter (in parentheses). In addition to the statistical errors, wc also show estimates laf the error budget from possible sources of systematic uncertainty. 
We consider a range of different systematic effects, which arc described in detail in § o; specifically wc evaluate the effect of the choice of dark matter 
halo model (ADM), adiabatic contraction (AAC), treatment of the background (ABackground) and the Solar wind charge exchange X-ray component 
(ASWCX), the extrapolation of the entropy model (AEntropy), dcprojection (A3d), the employed spectral model (ASpectral), removing the emissivity 
correction (AWeighting), varying the priors on the model parameters (AFit priors) the X-ray detectors used (AInstrument), the modelling of the stellar 
light (AStars), distance uncertainties (ADistance), fit radius (AFit radius) and covariance between the temperature and density data-points (ACovariance). 
We list the change in the marginalized value of each parameter for every test and, in parentheses, the statistical uncertainty on the parameter determined 
from the test. Note that the systematic error estimates should not in general be added in quadrature with the statistical error. 




a 

0) 



100 



800 
Radius (kpc) 



500 




Fig. 9. — Baryon fraction profile inferred from our best-fitting 
models (black line). The shaded grey region indicates the l-a sta- 
tistical uncertainty in the fits. The red shaded region indicates 
the gas fraction profile and l-a uncertainty. The dotted line indi- 
cates the best-ptiing Pnsmnlogiral value of fj, , based On the 5-year 
WMAP data (IDunkley et alJp009^), which lies significantly above 
the measured tj, lor most 01 tne radial range. We indicate the 
physical scales corresponding to Rvir and various other standard 
radii. 

of radius, and we tabulate the gas fraction [fn) and fb, 
measured at different over-densities, in Table E. Outside 
~100 kpc (~R25oo)! our model predicts that the fb pro- 



FlG. 10. — Marginalized probability density for fb.vir for different 
choices of the extrapolated entropy profile. The solid (black) line 
is for the default case, where we-extrapolate the measured entropy 
profile to Rvir. As shown in § pA, this extrapolation is consistent 
with trends seen in massive galaxy clusters. Also shown are the 
cases where we allow the entropy profile to break at '~80 kpc, and 
the logarithmic slope (7) is much steeper (7 =2.5; dashed red line) 
or flatter (7 = 0.0; blue dotted line). For clarity, we smoothed 
the distribution functions with a fourth-order Savitzky-Golay filter, 
which spans ~25% of the measured fb range. Although the besi- 
fitting values depend to some degree on this choice (see also § 5.3), 
in all cases fb,vir is close to the Cosmological value (0.17). 
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TABLE 4 

Baryon fraction results and error budget 



Test 


/s,2500 


fg, 500 


fg,vir 


/6,2500 


/i),500 


fb ,vir 


Marginalized 
Best-fit 


(0.024) 


"■""-0.02 

(0.053) 


^■^^-0.03 

(0.119) 


0.10 ±0.01 
(0.096) 




'^•^^-0.02 

(0.102) 


0.16 ±0.04 
(0.154) 



Systematic error estimates 



ADM profile 
AAC 

ABackground 
ASWCX 
AEntropy 
A3d 

ASpectral 
AWeighting 
AFit priors 
AInstrument 
AStars 

ADistance 

AFit radius 

ACovariance 



-0.004 (±0.003) 
-0.001 (±0.004) 
-0.003 (±0.003) 
±0 (±0.002) 
±0 (±0.004) 
+0.005 [t'ofot] 
(±0.003) 
-0.001 (±0.002) 
iHSI (±0.004) 
-0.008 (±0.005) 
±0 (±0.004) 

-0.002 
-0.007 
-0.002 



-0.024 (±0.01) 
-0.010 

-0.001 (±0.01) 

+0.001 (iHi) 
-0.018 (±0.01) 
+0.010 (±0.02) 

-0.026 (±0.01) 
±0 (±0.01) 

io;o2 (±0.02) 

-0.025 (±0.02) 
-0.006 (torn 



-0.085 (iHi) 
-0.019 



f+0.05 
\ -0.03 



+0.004 ' 
-0.003 , 
+0.007 ' 
-0.005 , 
'+O.OO4' 
-0.003 , 



+0.03 (±0.04) 

+0.03 {tVol) 
(±0.05) 
+0.04 (tH^s) 

to.o? (±0.05) 
+0.01 (±0.05) 

+0.02 f+0.06\ 
-0.04 \^-0.05 J 

-0.067 (±0.06) 
+0.008 (±0.05) 



-0.007 (tool 
-0.019 (±0.02) 
-0.005 (±0.01) 



rO.OOl 
-0.064 



fO.05 
-0.04 
f0.07 
-0.04 



+0.02 {tlT,] 



(±0.01) 
(±0.01) 

{-0.01) 

(±0.01) 

(±0.01) 

(±0.01) 

/ +o.oi\ 
V-0.01; 

(±0.01) 
(±0.01) 
(±0.01) 

/+o.o3^ 

1^-0.05 J 

(±0.01) 
-0.017 (±0.02) 
+0.002 (±0.01) 



-0.013 
-0.026 
-0.005 
+0.003 
±0 
+0.01 

+0.00 
-0.03 

-0.006 

+0.01 
-0.01 
-0.021 
+0.02 
-0.01 
+0.00 
-0.00 



-0.044 
-0.023 
±0 

+0.005 

+0.01 
-0.01 

+0.01 

+0.01 
-0.04 

-0.003 
±0.02 

-0.032 
-0.008 



(±0.01) 
(±0.02) 
(±0.02) 
(±0.01) 
(±0.02) 
(±0.03) 

(±0.02) 
(±0.01) 
(±0.02) 
(±0.02) 

'+0.03^ 
-0.02 I 
'+O.O2' 
-0.02 , 

'+O.O3' 
-0.02 , 



-0.006 
-0.032 
±0.00 (±0.02) 



-0.114 

-0.035 

+0.02 

+0.03 
+0.01 

-0.05 

+0.03 

+0.03 
-0.09 

f0.008 

+0.01 

-0.06 

-0.081 

+0.01 
-0.01 

-0.006 
-0.076 
f0.003 



(±0.01) 
(±0.04) 

-0.04\ 
-0.03 J 

-0.04 J 

(±0.05) 

f+0.06\ 
^^-0.08 J 

(±0.05) 
(±0.05) 
(±0.06) 
(±0.05) 

?+0.05" 
1^-0.06 

(±0.05) 

f+0.Q8\ 
\-0.0i J 

(±0.05) 



Note. — Marginalized values and 1-cr confidence regions for the gas fraction (/g,A) and baryon fraction {fb,^) measured at various 
ovcrdensitics (A). We also ^provide the best-fitting parameters in parentheses, and a breakdown of possible sources of systematic 
uncertainty, following Table pi We find that is reasonably robust to most sources of systematic uncertainty, especially within R2500. 



file rises gently, from '^10-15%. In Fig we show the 
posterior probability density function for fb^vir (fb extrap- 
olated to Rvir)j which is tightly constrained to around 
^15 %. very close to the Cosmological value. A poten- 



and supported by ob servations of galaxy clusters (e.g. 



Vikhlinin et al. 2006^, we also experimented with a so- 



called "cored logarithmic" mass model, whic h has been 



tiallj important source of uncertainty in this extrapola- 



tion is the shape of the entropy profile. In Fig p.O|we show 
the equivalent posterior probability density functions for 
two different entropy extrapolations, which should bound 
all l ikely entropy profiles. We discuss these fits in detail 



widely used in studies of elliptical galaxies (Binney & 
Trcmainc |20qj), including s ome recent stellar dynamical 



studies 
^0091) . 



leg 
This 



5.3. Clearly, uncertainties the entropy extrapolation 
do not produce a large systematic shift in our inferred 

fb, vir ■ 

5. SYSTEMATIC ERROR BUDGET 

In this section, we address the sensitivity of our re- 
sults to various data analysis choices that were made, 
including the choice of prior. In most cases, it is difficult 
or impractical to express these assumptions through a 
single additional model parameter over which one might 
hope to marginalize, and so we adopted the pragmatic 
approach of investigating how our results changed if the 
assumptions were arbitrarily adjusted. We focused on 
those systematic effects likely to have the greatest im- 
pact on our conclusions, and list in Tables |3| and ^ the 
change to the marginalized value of each key parame- 
ter for each test. We outline below how each test was 

perfqnUBd Thnsp rpa.Hprs unintprpstpH in ihp tnrhniral 



Thomas et al.||2007b| ; Gebhardt fc Thomas 
model tends to predict higher masses at 
large radii than NFW, resulting in a significantly larger 
Mvir, and correspondingly a smaller value of fb when 
extrapolated within Rvir- Nevertheless, at radii corre- 
sponding to higher mass overdensities, the extrapolation 
range is smaller and the discrepancy is significantly re- 
duced, so that the choice of DM model has little effect 
on fb,25oo ("ADM profile" in Tables | and |). In any 
case, we note that the cored logarithmic model is less 
well-motivated theoretically than the NFW profile. Fur- 
thermore, we find that the data slightly favour the NFW 
model; although the best-fitting is comparable be- 
tween the two cases, the ratio of the Bayesian evidence is 
~ 1.5 X 10~^, implying that the cored logarithmic model, 
with the adopted priors (a fiat prior on the asymptotic 
circular velocity, between 10 and 2000 km s~^, and a 
flat prior on logioVc, where Vc is the core radius, over the 
range < logiQrc < 3.) is a poorer description of the 
data at ^3.8-a. 
One modification to the (NFW) DM profile shape that 



detai lls of our analysis may wish to proceed directly to 
Sect ion |6| 

5.1. Dark Matter profile 

One of the major sources of uncertainty on the re- 
covered fb is the choice of DM mass model. While 
our default model (NFW) is theoretically motivated. 



is theoretically pred i cted is "adiabatic co ntraction" (Blu 
menthal et~al| |l986| ; [Gnedin et aI||2004D , where the DM 



halo density profile reacts to the gravitational influence 
of baryons that are condensing into stars by becoming 
cuspier. There is little observational evidence of this ef- 



fect (e.g. IH06I: Humphrey et al. 2009al : pnedin et al.||20G7| ; 

although, see Napolitano et al.| 2010|). and it may be over- 
estimated by existing models ( Abadi et al. 2009), so we 
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did not employ it in our default analysis. Nevertheless, 
we experimented with modifying the NFW halo profile 
to account for th i s effe ct with the algorithm described 
by iGncdin et all ( |2004D p|. The principal effect of this 
modification was to lower slightly the stellar M/L ratio 
and, consequently, the inferred fb ("AAC" in Tables |^ 
and|). 

5.2. Background 

Our treatment of the background is a potentially se- 
rious source of systematic error. To investigate this, we 
adoptrd an alternative approach in which we used, for 



the Chandra data, the standard blank-field events files 
distributed with the CALDB to extract a background 
spectrum for each annulus. Since the blank-field files 
for each CCD have different exposures, spectra were ac- 
cumulated for each CCD individually, scaled to a com- 
mon exposure time and then added. The spectra were 
renormalized to match the observed count-rate in the 9- 
12 keV band. These "template" spectra were then used 
as a background in Xspec, and the background model 
components were omitted from our fit. We found that 
this approach gave a poorer fit to the data than the full 
modelling procedure used by default, and the overall den- 
sity and temperature profiles were noticeably affected at 
large radii. For the Suzaku data, there are no blank- 
fields files available, but we used the standard non X-ray 
background spectra generated by xisnxbgen tool to re- 
move the instrumental background, but maintained the 
sky background model components in our fit. Fitting the 
resulting Chandra and Suzaku temperature and density 
profiles with our mass modelling apparatus, our results 
were not substantially affected ( "ABackground" in Ta- 
bles | and I). 

5.2.1. Solar Wind Charge Exchange 

An additional background component can arise from 
the interaction of the Solar wind with interstellar mate- 
rial and the Earth's exosphere. This should manifest 
itself as a time-variable, soft component that can be 
modelled as a series of narrow Gaussian lines, the in- 
tensi ty of which correlate with the Solar wind activity 
(e.g. 



ing t: 



Snowdcn ct al. 2004| ). In the interests of minimiz- 
le complexity of 



the spectral models, we do not by 
default explicitly account for this so-called "Solar wind 
charge exchange" (SWCX) in our fits. To some degree 
this component is degenerate with the sky background, 
which we fitted in our modelling, and so this omission 
may not be a significant source of bias in our analysis. 
Nevertheless, it is important to confirm that it is not. 

To assess the contamination from the SWCX compo- 
nent, we first identified periods of low Solar wind proton 
fiux, as measured by the Solar Wind Electron Proton 
Alpha Monitor (SWEPAM) instrument aboard t he Ad 



mas et al. 


19981) 




assumed t 


he S"V 





Following Snowdcn ct al. ( p004 ), we 



X component is negligible tor a So 
lar wind proton flux level <3 x lO^cm"^ s~^. Selecting 
only times for which SWEPAM measured a proton fiux 

Using the CONTRA code publicly available from 
http: / /www. astro. lsa.umich.edu/~ognedin/contra/ 

Based on the publicly released data available from 
http://www.srl.caltech.edu/ACE/ASC/index.html 



below this limit, we were left with ~145 ks of Suzaku 
data. Still, since the proton flux was only mildly en- 
hanced during the (short) excluded periods, we found 
that the soft (0.5-1.0 keV) X-ray count-rate was only 
enhanced by a few percent in the entire (i.e. unfiltered) 
data as compared to the quiescent period. We found that 
there was very poor SWEPAM coverage of the periods 
during which the Chandra data were taken, and so we did 
not attempt to filter the Chandra data in the same way. 
Instead we modified the sky background component by 
adding a number of narrow Gaussian lines (fix ed at the 
energ ies reported for the SWCX component by Snowden 



et al.). 

The corrected Chandra and Suzaku spectra were fit- 
ted to obtain the temperature and density profiles, and 
folded through our mass modelling apparatus. We found 
that our results were not substantially affected by cor- 
recting for the SWCX component ("ASWCX" in Ta- 
bles |and I). 

5.3. Entropy profile extrapolation 

The extrapolation of the entropy profile is one of the 
prime sources of uncertainty in ou r de termination of fb at 
large scales. As we will show in § |6.4 our default profile 
is, in fact , consistent with t rends seen in massive galaxy 
clusters ( |Pratt et aLl|2010| ), giving us confidence in its 
use. Nevertheless, it is important to explore plausible al- 
ternative shapes. Motivated by observed entropy profiles 
in galaxy groups and by the monotonically rising entropy 
profiles required by hydrostatic equilibrium, we investi- 
gated two pathological extrapolations that should bound 
any plausible model. Specifically, we adopted a powerlaw 
shape (entropy oc R"') outside ~80 kpc, with 7 = (i.e. 
a flat entropy profile), which is a secure lower bound- 
ary assuming stablity against convection, or 7 = 2.5. 
This latter value was chosen arbitrarily, and is far larger 
than has been observed in real systems. This choice pri- 
marily affects the extrapolated temperature, rather than 
the density profile, and so fb does not change substan- 
tiallv with these choices ( " AEntropy" in Tables and 0; 
Fig0). 

5.4. Deprojection 

In the present work, we opted to fit the projected, 
rather than the deprojected data (as done, for exam- 
ple, in HOG). In general, fitting the projected data leads 
to smaller statistical error b ars, but potentially large r 
systematic uncertainties (e.g. Gastaldello et al. 2007^), 
and so it is important to investigate the likely magni- 
tude of such errors. To do this, we examined the effect 
on our results of spherically deprojecting the data. We 
achieved this by using the Xspec projct mode]|^ To ac- 
count for emission projected into the line of sight from 
regions outside the outermost annulus, we added an apec 
plasma model to each annulus, with abundance 0.3 (con- 
sistent with the outermost annulus) and the temperature 
and normalization determined from projecting onto the 
line of sight the best-fitting gas temperature and density 

In practice, it was more convenient to emulate the behaviour 
of projct by adding multiple "vapec" plasma models in each annn-. 

with tht' relative normalizations tied appropriately (e.g. Kriss 
pt al. 19831). This allowed data from the multiple Suzaku instru- 
ments to DC fitted simultaneously. 
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4.2, but considering the models 



models described in £ 
only outside ~70 kpc. 

Although the error-bars were increased when using 
the deprojected (rather than the projected) profiles 
( "ADeprojection" in Tables @ and | ) , the best- fitting re- 
sults were not significantly changea. In Fig ^, we show 
a series of mass "data-points" obtained from the depro- 
jected density and temperature profiles, usin g the "tra- 
dition al" mas s analys is method described in Humphrey 
et al| (|2009ai see also [Humphrey fc Buote|p01(i| ). These 



data agree very well with the best-fitting mass model 
found in our projected analysis. Similarly, the entropy 
profile derived directly from the deprojected data agrees 
well with the results of our projected analysis (Fig ||). 

5.5. Priors 

Since the choice of priors on the various parameters is 
arbitrary in our analysis, it is important to determine to 
what extent they could affect our conclusions. To do this, 
we replaced each arbitrary choice in turn with an alter- 
native, reasonable prior. Specifically, for each parameter 
describing the entropy profile, we switched from a flat 
prior on that parameter to a flat prior on its logarithm. 
We replaced the flat prior on the stellar M/L ratio with 
a Gaussian prior, the mean and cr of which were deter- 
mined from the s tellar population synthesis model results 
reported in H06. We used a flat prior on the DM halo 
mass, rather than on its logarithm, and, instead of the 
flat prior on log cdm, we adopted the distri buti on of c 
aroun d M fo und by either Buote et al. ( 2007| ) or Maccio 
^t al. (2008) as a (Gaussian) prior. The effect of these 
choices is no larger than the statistical errors on each pa- 
rameter, especially for the baryon fraction measured at 
i?200 or higher overdensities ( " AFit priors" in Tables ^ 
and|). 

5.6. Stellar light 



As discussed in H06, careful treatment of the stellar 
light is essential for obtaining an accurate measurement 
of the gravitating mass proflle. Since the stellar light 
deprojection is not unique, in particular as the inclina- 
tion angle is not a priori known for an elliptical galaxy, 
we investigated the associated systematic uncertainty by 
first varying it within reasonable limits (§ |). Specifi- 
cally, we lowered it as far as 70°, since the highly ellip- 
tical projected isophotes of NGC 720 make it unlikely to 
be observed at a much lower inclination. As a more ex- 
treme test of the influence of the stellar light on our flts, 
we also experimented with excluding all the data within 
the central ~4 kpc, where the stars dominate the mass 
proflle. We flnd that the stellar M/L ratio becomes (un- 
surprisingly) much more uncertain when we adopt this 
approach, but the best-fitting masses and fb were not 
substantially affected ( "AStars" in Tables ^ and ||) . 

5.7. Emissivity correction 

In our default analysis, the projected temperature and 
density profile are weighted by the gas emissivity, folded 
through th e instrumental responses (for details, see Ap- 
pendix B of iGastaldello et alj^007l| ). Since the computa- 
tion of the gas emissivity assumes that the three dimen- 
sional gas abundance profile is identical to the projected 
profile (which is unlikely to be true), it is important to 



assess how sensitive our conclusions are to the emissivity 
correction. To do this, we adopted the extreme approach 
of ignoring the spatial variation of the gas emissivity al- 
together. We found that this had a very small effect on 
our results (AEmissivity in Tables ^ and |4|). 

5.8. Remaining tests 

We here outline the remaining tests we carried out, 
as summarized in Tables ^ and 13. First of all, since the 
inter-calibration of the Suzaku XIs units may not be per- 
fect, we experimented with using only one of the units in 
the Suzaku analysis, and cycled through each choice. We 
found that the XISl (back-illuminated) unit appeared 
most inconsistent with the Chandra data, and so also 
investigated using only units 0, 2 and 3. We also consid- 
ered using only the Chandra data. While these choices 
had a significant effect on fb,vir, we found that fb at higher 
overdensities was more resilient ( " AInstrument" ) . 

To assess the impact of various spectral-fitting data 
analysis choices on our results, in turn we varied the 
neutral hydrogen column density by ±25%, performed 
the fit over a restricted energy range (0.7-7.0 keV, or 
0.5-4.0 keV), and replaced the APEC plasma model 
with a MEKAL model. The impact of these choices is 
comparable to the statistical errors on the parameters 
("ASpectral"). 

To examine the error associated with distance uncer- 
tainties, we v aried the dis t ance to NGG 720 by thel-cr 
error given in 
suits is minor 



Tonry ct al.| ( ^OOl] ). The effect 
[ "ADistance" ) . To examine the sensitivity 
of our fits to the radial range studied, we excluded the 
data outside ~40 kpc. This exclusion made the extrap- 
olation more uncertain, especially for fb.vir ("AFit ra- 
dius" ) . Finally, to examine the possible errors associated 
with our treatment of the covariance between the den- 
sity data-points, we investigated adopting a more com- 
plete treatment that considers the covariance between 
all the temperature and density data-points, as well as 
adopting the more standard (but incorrect) approach of 
ignoring the covariance altogether. The effect is not large 
("AGovariance" in Tables | and |; Fig |). 

6. DISCUSSION 

The combination of the new, high quality data for 
NGG 720 and our updated, Bayesian hydrostatic mod- 
elling procedure has greatly improved our constraints on 
its gravitating matter and hot ISM. We here discuss the 
implications of these new measurements at some length. 

6.1. Hydrostatic equilibrium 

The ability of our hydrostatic model to reproduce the 
temperature and density profiles of the gas in NGC 720 
strongly suggests that the gas is close to hydrostatic; de- 
spite highly nontrivial temperature and density profiles, 
a smooth, physical mass model and a monotonically ris- 
ing entropy profile (required for stability against convec- 
tion) were able to reproduce them well. If the hydrostatic 
approximation is seriously in error, this would require a 
remarkable conspiracy between the temperature and den- 
sity data points. The closeness of the system to hydro- 
static is unsurprisi ng given its remarkablv relaxed X-ray 
morphology (Fig ^; [Buote et"aL 2002 ; Buote fc Canizares 
1996 . 1994 ) and la ck of substantial radio emission ( [Fab- 
piano 'ct~al.| 1987[ ). Numerical simulations of structure 
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formation suggest that deviations from hydrostatic equi- 
hbrium in morphologicaUy relaxed-looking systems are 
not large, so that errors in the recovered r nass distri- 
bution should be no larger than ^25 % (e.g . |Tsai et al 



1994 puotc k Tsai|[l995|: [Nagai ct aLl^OOTj ; pPiffaretti fc 
Valdarnm] ^008| ; |lkhg7rtal.|tj009| )r^cEa level ol non- 



thermal support is comparable to the ^10-20% discrep- 
ancy found betw een X-ray and stellar dynamical mass 



measurements by Churazov et al 
that are manifestly more disturbei 



008) for two galaxies 
tEan NGC 720. 
In the case of NGC 720. there is additional, albeit cir- 



cum^ tantial, evidence which supports the hydrostatic ap- 
proximation. First, the best- fitting virial mass is in good 
agreement with the modest (although uncertain) veloc- 
ity dispersion of the dwarf companions about NGC 720, 
(117 ± 54 km s^ ^, corresponding to a virial mass of 
4 ± 4 X IO^^Mq: |B rough et aH|2006| ). In our studies 
of other systems, we have used the good correspondence 
between the measured stellar M/L ratio and that pre- 
dicted by single-burst stellar p o pulat ion synthesis (SSP) 
models that assume a Kroupa ( ^001 ) IMF, to infer that 
there is little nonthermal pressure ( [Humphrey et aT 
^009aD . In the case of NGC 720, the measured M/Lk 



ratio (0.54 ± 0.05) is actually higher than the predic- 
tions of_tlic_SSPjmddsJ^iji±ii^^ 



way ^o bring these results into agreement is if the young 



f 



(-3 Gyr: [Humphrey fc Buotcj |2006| ; [Rcmbold et al.||2 Oq 
stellar population that dominates the light only repre 
sents recent star-formation superimposed on an older, 
underlying stellar population with a higher M/Lk ratio. 
Such a picture is supp orted bv the mu lt iple- a ged stellar 
population inferred by plembold et al, ( 2005| ). We note 
that the measured M/Lk value could also be reconciled 
with the predictions of SSP models if a Salpeter IMF 
is adopted (for which M/Lk'^0.54; taken in the context 
of our previous work this would indicate a non-universal 
IMF). Ahernatively, if AG operates, making NGC 720 
unique amongst the systems we have studied in the X- 
ray, the measured M/Lk would be 0.39 ± 0.04, also rec- 
onciling the fit result and the SSP predictions. While 
it is evidently difficult to draw strong conclusions based 
on the M/Lk ratio, none of these arguments indicates 
that the measured M/Lk is obviously underestimated. 
This is important since most models predict that non- 
hvdr ostatic effects lead to an underestimate of the mass 
(e.g. Nagai et al, 2007| [Churazov et al, 2008; Fang et al 
20091) 



•"Flumphrey fc BuOtc ' ( pOlOj ) studied NGC 720, as part 
of a sample of galaxies, groups and clusters, finding that 
the total gravitating mass profile in the central region 



(< 3( | kpc) is close in shape to a singular isothermal 
sphere distribution. This shape is what is approximately 
inferred from lens ing plus stellar dynami cs studies of sim- 
ilar galaxies (e.g. Koopmans et al. [2009 ). Therefore, any 
deviations from hydrostatic equilibrium are unlikely to 
have substantially affected the overall shape of the grav- 
itating mass profile, at least over this region. Neverthe- 
less, it is plausible that a constant fraction of the pressure 
over this region comes from nonthermal effects, which 
would just reduce the overall normalization of the mea- 
sured ma,ss mode] Such a, scenario would imply tha.t the 



putting it in more tension with theory. Conversely, it 
is unlikely that deviations from hydrostatic equilibrium 
could lead to the scale radius also being underestimated, 
since that would require the nonthermal pressure sup- 
port to increase dramatically outside ^-^20 kpc. Since 
NGC 720 is an isolated system, it is difficult to imagine 
what processes could give rise to such an effect. 

One effect which could produce deviations from hy- 
drostatic equilibrium at the smallest scales is gas rota- 
tion induced by spin-up of subsonically inflo wing gas. 
Cons i dering the very relaxed galaxy NGC 4649 Brighenti 



et al, ( 2009| ) demonstrated that angular momentum con- 



servation (if there is appreciable rotation in the stars) 
flattens the X-ray ellipticity significantly in the centre of 
the system (where deviations from hydrostatic equilib- 
rium are also most signific ant). NGC 720 exh ibits suffi- 
cient major-axis rotation (Binney et al. 1990), but only 
a hint of a central ell i pticit v rise within the innermost 



1 kpc ( |Buote et aL| ^002] ). Since our results are not 



very sensiti ve t o the exclusion of data within the central 



^4 kpc (§ p.8[ ), it seems unlikely that this effect could 
have an impact on our results, even if it is operating in 
NGC 720. 

Our conclusions on the state of hydrostatic e quilib 



rium in NGC 72 differ significantly from those of Diehl 
fc Static: (£007). Those authors investigated the X-ray 
and optical isophotal ellipticities at a fixed, small scale 
in an heterogenous sample of early- type galaxies, con- 
cluding that the lack of a correlation between them (in 
confiict with would be expected if hydrostatic equilib- 
rium holds exactly) indicates that hydrostatic equilib- 
rium is not ubiquitous and, therefore, should never be 
assumed (even approximately) in mass analysis (and, in 
particular, for NGC 720). We consider their conclusions, 
however, to be extreme and misleading. Firstly, their 
method does not allow them to say anything at all about 
individual objects, and only indicates that hydrostatic 
equilibrium is unlikely to hold perfectly in a galaxy ran- 
domly drawn from their sample, provided it is chosen 
without any consideration of its morphology. Since their 
sample contained galaxies with such large-scale asym- 
metries that we would not recommend the routine ap- 
plication of hydr ostatic mas s methods (e.g. NGC 4636 : 
Jones et al.||2002| and M84: [Finoguenov fc Jones] |200l| ), 
and since they focused on the smallest scales (where X- 
ray point-source removal is most challenging and where 
AGN-driven disturbances are most serious), in contrast 
to the large scales that are most important for dark 
matter analysis, the implications for morphol ogicallv re- 
laxed system s such as NGC 720 or NGC 4649 ( [Humphrey 
et al. 2008) are unclear. Secondly, even if it does not 



hold perfectly, hydrostatic equilibrium can still be a use- 
ful approximati o n. ^V ith their hydrodynamical models. 



Brighenti ct al. (200S) reconciled the X-ray and optical 
ellipticity profiles of NGC 4649 by requiring modest gas 
motions, but the gas remained very close to hydrostatic 
(especially outside the innermost '^l kpc). Substantial 
non-hydrostatic gas motions are likely to manifest them- 
selves as sharp features in the surface brightness and 
temperature profiles of th e hot gas, for examp l e those 



true IRvi,- is underestimated, while the scale radius is ap- et al.||201C| ). StiU, even for M 87, jChurazov et al 



proximately correct; hence Cvir would be underestimated, 



seen in the core of M87 (Forman et al._ 2007 ;_ Million 



found that the gas was close to hydrostatic away from 
the regions of most significant disturbance. In the case 
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of NGC 720, the profiles are significantly smoother than 
in M87 (Fig while the h y drost atic isophotal analysis 
undertaken by Buote et al. (2002) found evidence for a 
dark matter halo ellipticity in excellent agreement with 
the predictions of cosmological models (implying the gas 
is not grossly out of equilibrium). In summar y, in our de- 
tailed analys i s of NGC 720 and othe r systems ( plumphrey 
et al| ^009a|: iBrighenti et al.| ^0091 ), we find no evidence 
supporting Diehl fc Statler| 's contention that hydrostatic 
equilibrium is a very poor approximation in morphologi- 
cally relaxed galaxies (which seems little more than "guilt 
by association" ) . 

6.2. Baryon fraction 
6.2.1. Robustness of ih measurement 

Based on our self-consistent hydrostatic mass model for 
NGC 720, and the fit to the density profile, we were able 
to obtain unprecedented fb constraints in a galaxy-scale 
system at large radius. Since the current X-ray data only 
reach ~0.7R25oo (with the field of view of Suzaku being 
the primary limiting factor), most of the constraints are 
based on an extrapolation. Nevertheless, this is done in 
a well-motivated, self-consistent manner, relying only on 
how the mass and entropy profiles behave at large ra- 
dius, and the validity of the hydrostatic approximation 
at large scales. Of these, the dark matter halo mass pro- 
file has by far the biggest influence on the inferred fb , as 
discussed in § Specifically, fb,vir changes from ^0.15 
when an NFW dark matter halo is adopted to ^0.04 for 
a "cored logarithmic" dark matter mass model, which 
mostly reflects the higher Mvir in the latter case. While 
this is a concern, the NFW model is much better moti- 
vated cosmologically, and supported observationally by 
high-quality X-ray and lensing studies of other (albeit 



more massiv e) systems f Lewis ct al. 2003 : Pointecouteau 
fit al.| ^005| ; [Vikhlinin et~al| [GiroSzi et al"] ^007| ). 



Still, even if the true DM mass model differs unexpect- 
edly from NFW, we flnd that fb measured within regions 
of larger overdensity (e.g. i?250o) is much less sensitive to 
this issue, reflecting the smaller range of extrapolation. 

While the extrapolation of the entropy proflle in our 
analysis is also a potential cause for concern, this choice 
does not significantly affect our conclusions (Fig |lO|) . In 
fact, even allowing for a pathological variation in the 
slope of the entropy profile at large radius (7 ranging 
from 0-2.5), we find that the best-fitting fb,vir varies by 
only r^l-a. This apparent lack of sensitivity to 7 arises 
because changing the entropy distribution in our hydro- 
static models tends to affect the implied temperature 
more strongly than the gas density. 

Another factor which may complicate the extrapolated 
fb values i s th e hydrostatic approximation. As we ar- 



gued in § 6.1, there is good evidence to suggest the 
gas is hydrostatic within ~75 kpc, but we cannot di- 
rectly assess the validity of this approximation at larger 
scales. Most mechanisms that can produce strong de- 
viations from the single-phase, hydrostatic approxima- 
tion in an isolated galaxy woul d actually provide non- 



thermal support to the gas (e.g. IZappacosta ct al 



Nagai et aL||2007| ; [Churazov et al.| p008| ; [Fang et al.|p009[ ). 
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I'his would produce a flatter gas density profile than the 
purely hydrostatic case, thereby increasing the actual 
baryon fraction over that inferred from our extrapola- 



tion. (An exception is if the gas is globally outflowing, 
in which case the gas is over-pressured, and the hydro- 
static model would over-estimate fb. However, such a 
global outflow would have to be triggered by feedback 
from the central galaxy, so it seems implausible for there 
to be an outflow outside ^75 kpc, but no evidence of one 
at smaller scales.) Numerical cosmological simulations of 
hot halos around galaxies do, indeed, predi ct that they 



are quasi-hydrostatic (e.g. Grain et al 



Our measurement of fb given in 'I'a 



2OIOD . 

le ^ does not in- 



clude all of the baryons within the system, since it ig- 
nores the dwarf companions. In practice, the total K- 
band luminosity of the detec ted dwarf galaxies i s only 



40% of that of NGG720 ([Brough et alj P006D , and 



of these galaxies 
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If the lu- 
minosity function of these galaxies is similar to that of 
Milky Way dwarfs, their combined lumi nosity will be 



domi nated by the brightest few members ( Tollerud et al 
2008), and so the baryons hosted by as yet undetected. 



fainter dwarfs will not signiflcantly affect this estimate. 
Additional baryons could plausibly exist in a hard-to- 
detect extended stellar envelope, analogous to the "intra- 
cluster light" seen in clusters and massive groups, but 
galaxy formation models suggest that, for a system of 
this mass, at most ^^10% o f the total stellar m ass could 
be distributed in this way ( Purccll ct al. 2007). Finally, 
we note t hat NGG 720 itsel f does not contain s ignificant 
cool gas ( iHuchtmeiei] I1994 [Welch et alj|2010| ). Failing 
to include these uncounted sources of baryons, therefore, 
may bias our inferred fb values low by at most ^0.01. 

6.2.2. Implications 

Based on our hydrostatic mass models, the total stellar 
mass fraction within Rvh- is small (~ 0.035 ± 0.003), but 
in good agreement with that generally inferred f rom grav- 
itatioiial lensing studies of early -- type galaxies (Hoekstra 



t al. 2005; Gavazzi et al. 2007). Most of the baryons 
within Rvir are, however, in the form of the hot gas. In- 
cluding both gas and stars, the measured fb profile is rela- 
tively flat at ~0. 10-0. 15 between R2500 and Rvir (Fig^), 
becoming c ompletely consistent with the Gosmological 
value (0 .17: Dunklcy ct al. 2009) when extrapolated to 
Rvii-. In II06 we reported a lower value of fb within Rvir 
for NGC 720 (^0.04), based on simpler hydrostatic model 
flts to data from a much shallower Chandra observation. 
That constraint, however, was almost entirely driven by 
a strong prior that imposed a tight correlation between 
Mvir and fb; when relaxing this prior, we found that fb 
was very poorly constrained by those data. 

A number of authors have investigated the gas (or 
baryon) fraction in galaxy groups and clusters, typically 
flnding a shallow dependence on the virial mass of the 
halo. In the upper panel of Fig |ll| we show the gas frac- 
tion measurements within R2500 (/n.2R nn) versus M^snn 



for the group and clust er samples of Vikhlinin et al 
( |2006| ), iGastaldello et alj ( ^007b| ) and ISuirctaLl ( |2009|) F1. 



Overlaid is the locus of NGC 720, which is in good agree- 
ment of an extrapolation of the powerlaw trend inferred 

For clarity, we omit systems with large error-bars 
(o-(/j)> 0.01). Where th e samples overlap , we use the results, in 



order nf preference, from Gastaldello et al. (2007h) and Sun ct al 
(2009). 
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Fig. 11. — Top panel: Gas fraction within R2500 shown as 
function of M2500 for the isolated galaxy NGC 720 (fiTp.pnintprl 
star) and a literat 
Eooel^^iiaffionds 



ire sample of grniips and clusters ( Vikhlinin ct al 
Gastaldello et al. |2007bf six-pointed stars; pun 



et al' lpOOgj : triangles j. ine Dest powerlaw fit to the data is shown 
as a solid line. The gas fraction within R2500 of NGC 720 is clearly 
consistent with an extrapolation of the trend for galaxy groups 
and clusters. Centre panel: The same, but within an mrerdnnsity 
of 500. We overlay additional literature data from Dai et al. ( 2009 
circles). Lower panel: Baryon fraction within Rsijo as a function 
of M500 for NGC 720 aniJ the, mean data for the group and cluster 
sample of Gio^lini et al j 2Dp9t squares). The best-fitting powerlaw 
relation from 



Tiodim ct aT 



is overlaid (dashed line). NGC 720 is 
clearly consistent witti an extrapolation of this relationship to low 
masses. 

for the more massive systems. Fitting a model of the 
form 



logio fg = A + B\o^ 



10 



M 



2 X 10"Mp 



(3) 



to these data, using a method that takes into account er- 
rors on both axes and intrinsi c scatter in the y-directio n 
(specificany, that employed in Humphrey & Buote 2008), 
we found A = -1.08 ± 0.03 and B = 0.28 ± 0.04, with an 
intrinsic scatter of 0.11 dex. In the centre panel of Fig [ll|, 
we show a similar comparison between /g,5oo and M500, 
similarly finding that the locus of NGC 720 is consistent 
with an extrapolation of the trends for the massive sys- 
tems. Fitting a log- linear regression line (Eqn we 
obtained A = -1.02 ±0.02 and B = 0.15 ±0.03, with an 
intrin sic scatter of .06 d ex, in good agreement with a 
fit by piodini ct al\ ( p009| ) to a similar group and cluster 
sample, but omitting the NGC 720 datapoint. Since the 
temperature varies only by ^^1 order of magnitude for the 
almost 3 orders of magnitude mass range between the 



galaxies and clusters, hydrostatic equilibrium generally 
requires that massive systems have a more rapidly de- 
clining gas density profile than less-massive objects, and 
hence a more centrally peaked gas mass distribution, as 
compared to the distribution of dark matter. Given this 
trend, it is not surprising that lower-mass objects have a 
smaller fraction of their total gas mass within R2500 1 as 
indicated by the dependence of B on the overdensity. 

Also shown in F ig pi] (cent r e pari el) are a set of data- 
points taken from Dai et al. ( 2009| ), who fitted stacked 
Rosat All-Sky Survey (RASS) data of optically selected 
galaxy groups and clusters. These points appear quite 
discrepant with this relation (and the other data) be- 
low 3 X IO^'^Mq, where they fall sharply. This may 
indicate real differences in the properties of optically se- 
lected groups from those selected from X-ray surveys, 
but it may also reflect the significant systematic errors 
that are in volved with stacking the RASS data (see Dai 



et al. 2007) 



While the gas is the dominant baryonic component for 
the most massive clusters, the stellar-to-gas mass ratio is 
actually a strong function of Msoo- In NGC 720, approx- 
imately half of the baryons within R.»snn are in the stars. 



Accounting for both gas and stars piodini et aL ( 2009| ) 
obtained average fb constraints within R500 tor a galaxy 
group and cluster sample. In the lower panel of Fig [l^, 
we show their data as a function of M500, along with me 
best powerlaw fit to their data, and the locus of NGC 720. 
It is immediately clear that NGC 720 is in good agree- 
ment with the extrapola tion of the (shall ow*) trend they 
observe to low masses. Dai et al. ( 2009| ) also reported 
total baryon fractions for their galaxy group and cluster 
sample, finding rather lower values of fb.200 in their lowest 
mass systems (fb,2oo~0.04 at M2oo~ 3 x lO^^Mg). This 
is significantly below fb,200 observed in NGC 720 (which 
has M200 an order of magnitude lower) , but is compara- 
ble to its stellar mass fraction. Still, as discussed above, 
the gas fraction estimates from their work appear dis- 
crepant with studies of X-ray selected groups and this 
may explain the inconsistency. 

Our study of NGC 720 clearly demonstrates that it is 
possible for a system with Mvir as low as a few times 
1O^^M0 to retain a massive, quasi-hydrostatic hot halo. 
The total baryon fraction, when extrapolated to the 
virial radius, is consistent with the Cosmological value, 
indicating that massive dark matter halos around galax- 
ies could be a significant re servoir for some fraction o f 
the "missing baryons" (e.g. [Fukugita fc Pcebles| |200^ ) . 
The biggest spiral galaxies are believed to reside m ha- 
los of comparable mass, making it particularly interest- 
ing to compare the hot halo around NGC 720 to those 
predicted around spiral galaxies. Intriguingly, the total 
gas content of the system (gas fraction within R200 of 
0.08 ± 0.03) does agree well with recent nu merical disk 



galaxy simulations at this mass range (e.g. Crain et al 
2010). In these simulations, most of the gas is in the 



hot phase, suggesting that the cool baryon content of 
the most massive spiral galaxies (i.e. considering only 
stars and cool gas) should be systematically lower than 
fb measured in NGC 720. 



McGaugh et al. ( |2010| , see also [Dai et al.| |2009| ) has 

assembled fb.500 estimates for a sample of massive disk 
galaxies, omitting any putative hot halo from the esti- 
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mate. Indeed, fb.500 in these systems does lie slightly 
below that in NGC 720. Still, the steUar mass fraction 
in these spiral galaxies is substantially higher than in 
NGC 720, and the modest discrepancy in fb,5oo leaves 
little room for a massive hot halo around the spirals. 
This appears inconsistent with the simulations and, at 
face value, might suggest very different star formation 
efficiencies between the different types of galaxy. One 
difficulty with this interpretation, however, is the young 
('--^3 Gyr) dominant stellar population in NGC 720, sug- 
gesting a recent formation epoch, possibl y from a merger 



event involving massive disk galaxies ( Rembold et al 
p005| ). Since the gas fraction can only be reduced during 
a merger (as it may be ejected from the galaxy, funnelled 
onto the central black hole or converted into stars), this 
implies there must have been at least as much gas as 
stars within R500 in the progenitors. Even allowing for 
stellar mass loss over the last ~3 Gyr, this would require 
similarly inefficient star formation in the precursor (spi- 
ral?) galaxies. \ Ve caution that t he M500 values for the 
spiral galaxies in McGaugh et al. were obtained by scal- 
ing the rotation velocity, rather than formal modelling, 
which could plausibly introduce some uncertainties into 
both M500 and fb.500 reported in that work. 

Finally, it is interesting to compare our results to mea- 
surements of fb in the Milky Way, around which there is 
indirect evidence for an extended, hot halo. This system 
has M^w in the range ~ 1-3 x lO^^M© (|Klypin et aL[|2002 



pakamoto et al.||2003[ ), close to that of NGC 720 . Base d 



on the baryonic mass estimate of |Flynn et al.| (2006), 
which does not consider any putative hot halo, ib,vir hes 
in the range ~0. 02-0. 06. This is suggestively below our 
measurement for NGC 720 and the discrepancy is most 
acute towards the more massive end of the Milky Way 
mass range (i.e. almost the same mass as NGC 720). 

6.3. Dark Matter Halo 

We have obtained a high-quality measurement of the 
gravitating mass profile for the isolated elliptical galaxy 
NGC 720 between -0.002Rvir and 0.2Rvh. (-0.7R25oo)- 



ies fc.g. Buote et al. 12002: Humphrcv ct al 


2006 


. 2009a: 


Fukazawa ct alj 2006|: P'SuUivan et al. 


to07 


Zhang et al. 


p007: Nagino & Matsushita 2009), the data are inconsis- 



quality constraints on the mass profile in NGC 720, this 
means that dark matter is required at 19.6-(t. 

With these data from deep Chandra and Suzaku ob- 
servations, we have been able to obtain unprecedented 
constraints on the virial mass and dark matter halo con- 
centration at this mass range (Mvir= 3.1±0.4x lO^^M©). 
Such a low mass (only a few times more massive than the 
Milky Way) is consistent with our picture of this sys- 
tem as an isolated elliptical galaxy, rather than a galaxy 
group per se. These results are in agr eement with, but 



very m uch tighter than, those found by Humphrey et al 



(2006) for the same system. This improvement partially 
reflects the much better data used here (^100 ks of Chan- 
dra data, plus ^180 ks of Suzaku ti me, as compared t o 
the ~17 ks Chandra dataset used by Humphrey et al, ), 
but " ' - . ^ 



m aiao 011 account of varlouH I'cfiiieineiiLy made Lo our 
mass modelling procedure in the interim. Of particular 
significance is the adoption of parametrized models for 
the entropy (rather than the temperature) distribution. 



allowing us to enforce a monotonically rising entropy pro- 
file (i.e. the Schwarzschild criterion for stability against 
convection). This efficiently eliminates unphysical por- 
tions of parameter space. Similarly, explicitly requiring 
the gas to remain approximately hydrostatic out to the 
virial radius also restricts parameter space to regions con- 
taining only physical models, whereas the fitting proce- 
dure of H06| was less restrictive. 

In large part, the tight constraints on the virial mass 
and concentration are due to the ability of the data 
to constrain the flattening of the mass profile outside 
~20 kpc. This scale is much larger than the effective 
radius of the galaxy (3.1 kpc), and is thus unrelated to 
the s tellar-dark matter conspiracy ( Humphrey fc Buote 
2010), and instead must indicate the scale radius of the 
dark matter halo. As discussed in Gastaldcllo et al 



(2007b), the ability to constrain the virial quantities is 
very sensitive to the accurate measurement of the scale 
radius. This flattening is not model dependent; when 
fitting the "cored logarithmic" DM model, we similarly 
found a flattening in the total gravitating mass proflle at 
this scaleFl. 

The srnall DM scale radius in NGC 720 indicates a 
high concentration; indeed, we find that it lies 3 
times the intrinsic scatter above the predicted median 
relation found in dark matter only simul ations for the 
"ACDMl" ( "concordance" ) cosmology by Maccio et al 



relat 



ion 



(|2008|). [Buote ct alj (|2007|) measured the c-M 
empirically tor a sample of early- type galaxies, groups 
and clusters, finding hig her concentrat ions in the galaxy 
regime (< 10"Mq) than [Maccio et aL| StiU, we find that 
NGC 720 lies ~2 times the intrinsic scatter above this 
relation^. The extreme isolation of NGC 720 suggests 
an early epoch of formation (such that the brightest sub- 
halos have merged to form the central object), which 
should give rise to a higher DM halo concentration, and 
may provide a possible exp lanation for the measured Cvir 
(e.g. Zentner et al. 2005). The very young (^^3 Gyr) 
mean age for the central stellar population may simply 
reflect stars formed in a recent gas-rich merger rather 
than a realistic estimate of the age of the population as 
a whole (see discussion in § 6.1). 

One factor which could lead to an over-estimate of Cvir 
i s a failure to model correctly the stellar mass com ponent 



(Mamon & Loka; 
young stellar popu 



2005; Humphrey ct al. 2006). If the 



ation which dominates the light is, in- 



deed, only a fraction of the total population in the galaxy, 
it is plausible that the assumption that mass follows light 
is not correct (although the modest colour gradient of the 
galaxy is not consistent with a dr amatic violation of this 
assumption: [Peletier et al. 1990). However, excluding 



We nntp thot tho 1-.f>gt-fittiriir magg profile agrees with that 
reported by Humphrey & Buote ( |2010[ ), who found that, within 
~30 kpc, it IS approximately consistent with a powerlaw. In that 
case, the data lacked the radial coverage to constrain the subtle 
flattening between '~20-30 kpc, within which there is only a single 
data bin. Furthermore, the slope of the mass proflle was insensitive 
n thn oyact radial fitting range adopted, SO our conclusions in 
Humphrey & Buotn (GOlO) are not affected by our failure to correct 



tor tnis nattcmnL 

Strictly speaking, this comparison is not independent ,_aiiic^ 
Nnr^ 790 i»raQ one of the systems used to determine the 
et al 



less, 



Buote 



(pOOTj) relation, albeit with much poorer data. Nevertne- 
•JGu i20 does not appear to be unusual in comparison to the 
other systems used by these authors, and so it is unlikely that this 
relation is largely driven by this one data-point. 
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As expected for gas which is approximately hydro- 
static, we obtain a good fit to the temperature and den- 
sity profiles with a monotonically rising entropy profile. 
Although the entropy profile slope in the inner part of the 
system (/3i) is close to the canonical f .1 (Table 0) pre- 
dicted by gravitational structure formation simulations 
(considcririg: on ly gravity without cooling or feedback: 
Voit ct al. [2005 ) , the profile flattens at smaller and larger 
radii. Thi s is consistent w i th two of the three galaxies 
studied in Humphrey et al, ( 2009a ), and similarly shaped 
profiles have also been observed in aalaxv groups fe.g. 
Jetha et al.| ^ 007 ; Finoguenov et al. 2007; ICastaldcllo 



"alM 

ETFT 



OOTi: iSiiTet al.||2009| [Flohic et aLfpOfOU 

T% we show the entropy prohle as a function 
and scaled by the "characteristic entropy" , K 



et 



Fig. 12. — 1-(T confidence region for the entropy profile model 
of NGC 720, scaled by its characteristic entropy K500 (=27.9 keV 
cm^), and shown as a function of Rmq. We overlay the deprojected 
entropy data-points (shown in Fig H), similarly scaled. The dotted 
line ^nrHr-ntoQ tho "hngf^linr." prorHr-tinn from gravitational structurc 
formation (Voit et al. 2005). Also shown is the scaled entropy 
model, corrected tor tne gas fraction profile ("/j corrected"; see 
text), which agrees very well with the baseline model. 




R (kpc) 

Fig. 13. — Cooling time of the gas (solid line), and its 1-cr confi- 
dence region (shaded region). Error bars do not incorporate abun- 
dance errors (which are less than ~20%). Gas inside ~60 kpc has 
a cooling time shorter than the Hubble Time (shown as a dashed 
line). 



the central ^4 kpc ('^l.SRp) does not significantly af 



feet the concentration (§ 5.8). Alternatively, if adiabatic 



contraction operates, it should increase the cuspiness of 
the central DM halo profile, for a fixed stellar M/L ra- 
tio, which could also lead to c „ir being over-est i mate d 
if it is not accounted for (e.g. [Napolitano et al" |oTo). 
Fitting an AC model (§ 5.1) also led to a slightly less 
concentrated halo, but the effect was not large enough 



t o rec oncile the best-fitting value with the Maccio et al 
( 2008 ) predictions. This comparative lack of sensitivity 



of the concentration to the details of the stellar mod- 
elling is not altogether surprising, since Re (3.1 kpc) is 
far smaller than the scale radius of the DM halo. 

6.4. Entropy profile 



of Rj 



(,e-g. 



500 



. Pratt ct al. 2010 ). Clearly the profile is significantly 
flatter than the "baseline" model predicted by gravity- 



only simulations (Voit et al. [2005 ), and the normaliza- 
tion is significantly enhanced, reflecting the injection of 
entropy by non-gravitational processes, especially in the 
inner regions . In comparison with galaxy gro ups and 
clusters (e.g. iPratt et al.|p010| ; |Sun et al.|^009| ), the en- 
tropy proflle IS more discrepant with the baseline model, 
reflecting the greater impact of nongravitational heating 
on a lower mass halo. Nevertheless, the entropy injection 
still appears insufficient to evacuate a significant fraction 
of the baryons from the halo, as indicated by the approxi- 
mate baryonic closure within Rvir- While this may reflect 
a finely-balanced heating mechanism in the system, we 
note that most of the gas lies at large radii, where the 
entropy profiles are much less discrepant (for example, 
^65% of the gas lies outside Rsoo)- 

If nongravitational processes primarily redistribute the 
gas, rather than raising its temperature, we might expect 
the entropy profile to be simply relat ed to the baseline 
model and the gas mass profile. Indeed Pratt et al. (2010) 
demonstrated with their cluster sample that the prod- 
uct S{R) {fg{R)/fb,u)^^^, where S is the entropy profile, 
fg{R) is the gas mass fraction profile (Fig ^ and fb,u is 
the Universal baryon fraction (0.17), is very close to the 
baseline model. It is interesting to investigate whether 
this also holds for lower-mass halos, where the impact 
of nongravitational heating is systematically larger. We 
show the "/g corrected" entropy profile for NGC 720 in 
Fig ^ finding that it also lies very close the gravity-only 
predictions, even in a ^Milky Way-mass halo. We note 
that the corrected entropy profile agrees with this rela- 
tion even when extrapolated outside the projected radii 
at which there are data, which provides support for our 
extrapolated entropy profile. 

In order to interpret the impact of non-gravitational 
processes on the entropy distribution, it is helpful to con- 
sider the cooling time of the gas. In Fig |l^, we show 
that, within ~60 kpc, it is shorter than the Hubble time. 
Nevertheless, only ~6% of the gas is currently within 
this region, implying that most of the gas in NGC 720 
will not have had sufficient time to cool. In part this re- 
flects the significant non-gravitational heating that has 
raised the entropy even in these outer regions signifi- 
cantly above the baseline model (Fig |l3) . The average 
stellar age in NGC 720 is young (~3 Gyr), and if this re- 
flects the time since NGC 720 formed, only the gas within 
~15 kpc would have had time to cool, within which there 
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is presently only ^^0.5% IO^Mq) of the total gas mass. 

6.5. Abundance gradient 

One of the intriguing results from our study is the 
detection, with both Chandra and Suzaku, of a strong 
negative abundance gradient in the hot gas of NGC 720 
(Fig H), making it arguably the lowest-mass system in 
which such a gradient has been definitively detected. Al- 
though the metallicity of the ISM in early-type galax- 
ies with intermediate to high X-ray luminosities is now 
known to be ^Solar, and comparable to the stella r metal- 



licity (Humphrey & Buote 2006; Ji et al. ^009 ), simple 



chemical enrichment models typically predict substan- 
tially more metals in the central part of the halo, where 
enrichment fro m supcrnovae and ste llar mass-loss is con- 
centrated (e.g. Pipino et al. 2005). Solutions to this 
problem may involve large-scale mixing due to AGN-ISM 
inte raction, or the a ssemb ly of s truct ure through merg- 
ers ([Mathews et al.| |2004t |Cora| |2006| ; |Cox et alj ^006|) . 
Clearly the spatial distribution of the metals in the hot 
gas is a useful constraint on any model predicting large- 
scale mixing. 

While, to date, very little is known about the radial 
dependence of Zpo in the ISM of galaxy-scale halos, re- 
cent Chandra and XMM studies of relaxed galaxy groups 
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were able to explain the overall shape of these group- 
scale profiles by invoking large-scale flows driven by low- 
level AGN heating. It remains to be seen whether such 
a model could also explain the observed abundance pro- 
file of NGC 720 which, residing in a much shallower po- 
tential well, should be more susceptible to the effects of 
feedback. 

7. CONCLUSIONS 

Based on an analysis of data from deep new Chan- 
dra and Suzaku observations, we have obtained un- 
precedented constraints on the temperature, density and 
heavy metal abundance in the hot gas halo of the Milky 
Way-mass isolated elliptical galaxy NGC 720. This has 
allowed us to measure in fine detail both the total grav- 
itating mass distribution and the baryon fraction out to 
large scales. In summary: 

1. The relaxed X-ray morphology of NGC 720, lack of 
radio emission, and the good fits obtained from our hy- 
drostatic models to the nontrivial temperature and den- 
sity profile indicate that the gas must be close to hydro- 
static. 

2. A constant M/L ratio model for the gravitating 
mass distribution is ruled out in favour of a model in- 
cluding an NFW dark matter halo at high significance 
(-20-cr). 

3. Assuming an NFW dark matter halo, we ob- 
tained unprecedented, tight constraints on both the virial 
mass and concentration for such a low-mass system 



(M25oo=l-6 ± 0.2 X IO^^Mq with systematic errors typ- 
ically <20%; Mvir= i-ltoi X lO^^M©). This confirms 
that an isolated elliptical galaxy that is not at the cen- 
tre of a massive group can maintain a substantial dark 
matter halo. 

4. The total gas mass inferred within Rvir substan- 
tially exceeds the stellar mass, with most of the gas lying 
outside '-^lOO kpc. This supports theoretical predictions 
that Milky Way-mass galaxies can host massive, quasi- 
hydrostatic hot gas halos. It remains to be seen if less 
isolated galaxies can, similarly, maintain such a halo. 

5. Self-consistently extrapolated to Rvir, the baryon 
fraction is consistent with the Cosmological value (0.17), 
indicating that it is possible even for a ^Milky Way-mass 
galaxy to be baryonically closed, at least if (relatively) 
isolated. 

6. Within both R2500 ^^'^ R-5001 the gas fraction (and 
the baryon fraction in the case of R500) of NGC 720 are 
both consistent with an extrapolation to low masses of 
the trends seen for groups and clusters. 

7 After correcting for the gas fraction, the entropy pro- 
file is close to the self-similar predictions of gravitational 
structure formation simulations, as observed in massive 
galaxy clusters. 

8. Both the Chandra and Suzaku data reveal evi- 
dence of a strong heavy metal abundance gradient, qual- 
itatively similar to those observed in relaxed, massive 
galaxy groups and clusters. 
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